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Abstract 

We provide an introduction to the basic concepts of chiral perturbation theory and discuss some 
recent developments in the manifestly Lorentz- invariant formulation of the one-nucleon sector. 
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1 Introduction 

Effective field theory (EFT) is a powerful tool in the description of the strong interactions at low 
energies. The central idea is due to Weinberg |Weinberg, 1979] : 

"... if one writes down the most general possible Lagrangian, including all terms consistent 
with assumed symmetry principles, and then calculates matrix elements with this Lagrangian 
to any given order of perturbation theory, the result will simply be the most general possible 
S-matrix consistent with analyticity, perturbative unitarity, cluster decomposition and the 
assumed symmetry principles." 

In general, an EFT is an approximation to a (more) fundamental theory, designed to be valid in a 
certain kinematical domain. Instead of solving the underlying theory, the processes under investigation 
are described in terms of a suitable set of effective degrees of freedom, dominating the phenomena 
in the particular energy region. In the context of the strong interactions, the underlying theory is 
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Fundamental theory 


Effective field theory 


Theoretical framework 


QCD 


ChPT 


Degrees of freedom 


Quarks and gluons 


Goldstone bosons (+ other hadrons) 


Parameters 


g 3 + quark masses 


Low-energy coupling constants + quark masses 



Table 1: Comparison of QCD and ChPT. 



quantum chromodynamics (QCD) — a gauge theory with color SU(3) as the gauge group. Under normal 
conditions, the fundamental degrees of freedom of QCD, namely, quarks and gluons, do not show up as 
free particles. One assumes that any asymptotically observed hadron must be in a color-singlet state, 
i.e., a physically observable state is invariant under SU(3) color transformations. The strong increase 
of the running coupling for large distances possibly provides a mechanism for the color confinement. 

For the low-energy properties of the strong interactions and the setting up of a corresponding EFT 
description, another phenomenon is of vital importance. The masses of the up and down quarks and, to 
a lesser extent, also of the strange quark are sufficiently small that the dynamics of QCD in the chiral 
limit, i.e., for massless quarks, is believed to resemble that of the "real" world. Although a rigorous 
mathematical proof is not yet available, there are good reasons to assume that a dynamical spontaneous 
symmetry breaking emerges from the chiral limit. Examples of indications for this to happen are the 
comparatively small masses of the pseudoscalar octet, the absence of a parity doubling in the low-energy 
spectrum of hadrons, and a non- vanishing scalar singlet quark condensate. 

According to the Goldstone theorem, a breakdown of the chiral SU(3) L x SU(3) fl symmetry at the 
Lagrangian level to the SU(3) y symmetry in the ground state implies the existence of eight massless 
pseudoscalar Goldstone bosons. The finite masses of the pseudoscalar octet of the real world are 
attributed to the explicit chiral symmetry breaking due to the quark masses in the QCD Lagrangian. 
Due to the vanishing of the Goldstone boson masses in the chiral limit in combination with their 
vanishing interactions in the zero-energy limit, a derivative and quark-mass expansion is the natural 
scenario for an EFT. The corresponding method is called (mesonic) chiral perturbation theory (ChPT) 
|Gasser and Leutwyler, 1984] , with the Goldstone bosons as the relevant effective degrees of freedom 
(see Tabled]). 

Using these effective degrees of freedom, physical quantities are calculated in terms of an expansion 
in q/A x , where q stands for momenta or masses of the pseudoscalar octet that are smaller than the 
energy/mass scale A x = 0(1 GeV) associated with spontaneous symmetry breaking. Since an EFT is 
based on the most general Lagrangian, which includes all terms that are compatible with the symmetries 
of the underlying theory, the corresponding Lagrangian contains an infinite number of terms, where each 
term is accompanied by a low-energy coupling constant (LEC). The method that allows one to decide 
which terms contribute in a calculation up to a certain accuracy is called Weinberg's power counting. In 
the mesonic sector, the combination of dimensional regularization with the modified minimal subtraction 
scheme of ChPT leads to a straightforward correspondence between the loop expansion and the chiral 
expansion in terms of momenta and quark masses at a fixed ratio. In actual calculations only a finite 
number of terms in the expansion in q/A x has to be considered and thus one has predictive power. 
What distinguishes the EFT approach from purely phenomenological approaches is the possibility of a 
systematic improvement. Mesonic ChPT has been tremendously successful and may be considered as 
a full-grown and mature area of low-energy particle physics. 

The situation gets more complicated once other hadronic degrees of freedom beyond the Goldstone 
bosons are considered. Together with such hadrons, another scale of the order of the chiral symmetry 
breaking scale A x enters the problem and the methods of the pure Goldstone-boson sector cannot be 
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transferred one to one. For example, in the extension to the one-nucleon sector the correspondence 
between the loop expansion and the chiral expansion, at first sight, seems to be lost: higher- loop 
diagrams can contribute to terms as low as 0(q 2 ) |Gasser et al., 1988] . For a long time this was 
interpreted as the absence of a systematic power counting in the relativistic formulation of ChPT. 
However, over the last decade new developments in devising a suitable renormalization scheme have 
led to a simple and consistent power counting for the renormalized diagrams of a manifestly Lorentz- 
invariant approach. 

The purpose of this article is to first provide a pedagogical introduction to the basic concepts 
of ChPT and to then present the more recent developments of a manifestly Lorentz-invariant ap- 
proach to the one-nucleon sector. It is definitely not intended to give a survey of the vast literature 
on ChPT and its various extensions in terms of chiral effective field theories. For further informa- 
tion the interested reader is referred to review articles and lecture notes addressing different topics 
with various priorities |Bijnens, 1993] , |Georgi, 1993j , Ecker, 1995 , Pich, 199~5] , Bernard et al., 1 995 1 , 
|Hemmert et al., 1998| , |Burgess, 2000 



|Bijnens, 20071 , |Bernard and Meifiner, 20071 , |Bernard, 20"08 



Scherer, 2003 , Scherer and Schindler, 2005] , |Epelbaum, 2006 



The article is organized as follows. Section [2] contains a discussion of the chiral symmetry of QCD, 
spontaneous symmetry breaking, and the Goldstone theorem. In Sec. |3J the basic concepts of mesonic 
ChPT are developed. Section [5] is devoted to baryonic ChPT. The power-counting problem is illustrated 
and solutions in terms of suitable renormalization conditions are presented. Section [5] contains a few 
selected applications of the manifestly Lorentz-invariant approach to nucleon properties. 



2 Chiral symmetry and spontaneous symmetry breaking 

The essential ingredients to setting up chiral perturbation theory as the effective field theory of the 
strong interactions are the chiral SU(3) L x 811(3)^ symmetry of QCD for massless u, d, and s quarks 
and the emergence of a spontaneous breakdown to the vectorial subgroup SU(3)y. 



2.1 Quantum chromo dynamics and chiral symmetry 

QCD is the gauge theory of the strong interactions |Gross and Wilczek, 1973| , |Weinberg, 19 73 1, 
|Fritzsch et al., 1973] with color SU(3) as the underlying gauge group. Historically, the color degree 
of freedom was introduced into the quark model to account for the Pauli principle in the description of 
baryons as three-quark states. The matter fields of QCD are the so-called quarks which are spin-1/2 
fermions, with six different flavors (u, d, s, c, b, t) in addition to their three possible colors (see Table 
[2]). Since quarks have not been observed as asymptotically free states, the meaning of quark masses 
and their numerical values are tightly connected with the method by which they are extracted from 
hadronic properties (see Ref. Manohar and Sachrajda, 2008 for a thorough discussion). 



2.1.1 The QCD Lagrangian 

The QCD Lagrangian can be obtained from the Lagrangian for free quarks by applying the gauge 
principle with respect to the group SU(3) of all unitary, unimodular, 3x3 matrices. Denoting the 
quark field components by qj,A,a, where / = 1, ■ • • , 6 refers to the flavor index, A = 1,2,3 to the color 
index, and a = 1, • • • ,4 to the Dirac spinor index, respectively, the "free" quark Lagrangian without 
interaction may be regarded as the sum of 6 x 3 = 18 free fermion Lagrangians: 

6 3 4 

Aree quarks = X X X ~ m /^aa')<?/,Aa'- (1) 

/=1 A=l a,a'=l 
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Flavor 


u 


d 


S 


Charge lei 

O L J 


2/3 


-1/3 


-1/3 


Mass [MeV] 


1.5-3.3 


3.5-6.0 


70 - 130 


Flavor 


c 


b 


t 


Charge lei 

o L J 


2/3 


-1/3 


2/3 


Mass [GeV] 


1.27+831 


a 9n+ - 17 


171.2 ±2.1 



Table 2: Quark flavors and their charges and masses. See Manohar and Sachrajda, 2008 for details. 



Suppressing the Dirac spinor index and introducing for each quark flavor / a color triplet 

/ 

Qf = 



QfA 

m h (2) 



V m 

the gauge principle is applied with respect to the group SU(3), i.e., all qf are subject to the same local 
SU(3) transformation: 

/ 8 A c \ 

q f ^q' f = eX p ( -i ^ ©ay) qf = Uq f , (3) 

where the eight denote Gell-Mann matrices acting in color space and the a are smooth, real functions 
in Minkowski space. Whenever convenient, we will make use of the summation convention implying a 
summation over repeated indices. Introducing eight gauge potentials A af j,, transforming as 

A, = .^3 A' = UA^ + -d,UU\ (4) 

the covariant derivative of the quark field, by construction, transforms as the quark field: 

Di*Qf = {.% + igsA^qf i (D^qf)' = D'^q' f = UD^q f . (5) 

In Eq. (JHJ), (73 denotes the strong coupling constant. In order to treat the gauge potentials as dynamical 
degrees of freedom, one defines a generalization of the field strength tensor to the non-Abelian case as 



rafiu 



d^Aau d u A a ^i Q^fabcAbfiAcvi (6) 



where, suppressing the superscript c in the Gell-Mann matrices, the standard totally antisymmetric 
SU(3) structure constants are given by (see Table [3]) 

/a bc = ^Tr([A a ,A b ]A c ). (7) 
At 

Given Eq. (J3J), the field strength tensor transforms under SU(3) as 

£v = Go^y >-> UG^uK (8) 

The QCD Lagrangian obtained by applying the gauge principle to the free Lagrangian of Eq. (CQ), finally, 

reads 

£qcd= E q f ^iP-m f )q f --g a , u g^. (9) 

r U,d,S, 

J~ c,b,t 



abc 


123 


147 


156 


246 


257 


345 


367 


458 


678 


fabc 


1 


i 

2 


i 

2 


i 

2 


i 

2 


i 

2 


i 

2 




|V3 



Table 3: Totally antisymmetric non-vanishing structure constants of SU(3): [4f, ^] = ifabc^f- 



From the point of view of gauge invariance the strong-interaction Lagrangian could also involve a 
term of the type 

£e = ^W^^r> e oi23 = 1, (10) 

where e^pu denotes the totally antisymmetric Levi-Civita tensor. The so-called 8 term of Eq. (TTUT) 
implies an explicit P and CP violation of the strong interactions which, for example, would give rise to 
an electric dipole moment of the neutron. The present empirical information indicates that the 9 term 
is small and, in the following, we will omit Eq. (flOl) from our discussion. 

2.1.2 Chiral limit 

The terminology chiral limit refers to massless quarks, resulting in an important additional global 
symmetry of the QCD Lagrangian which will be discussed in the following. We introduce the chirality 
matrix 75 = 7 5 = 27°7 1 7 2 7 3 = 75, {7^,75} = 0, 75 = 1, and define the projection operators 

P L =\{l-ls) = Pl P« = i(l + 75 )=4- (11) 

These operators satisfy the completeness relation Pl + Pr = 1, are idempotent, P| = Pl, P| = Pr, 
and respect the orthogonality relations PlPr = PrPl = 0. When applied to the solutions of the 
free massless Dirac equation, the operators Pr and Pl project to the positive and negative helicity 
eigenstates, hence the subscripts R and L for right-handed and left-handed, respectively. 
Omitting color and flavor indices, we introduce left- and right-handed quark fields as 

<?l = Pm and QR = PrQ- (12) 

A quadratic form containing any of the 16 independent 4x4 matrices {1, 7^, 75, 7^75, a^} can be 
decomposed as 

- r f ViFiiL + Qr^iQr for Ti G {YiYls} m \ 
q lQ { qRT 2 qL + qLT 2 qR for r 2 G{l, 7 5,^} ' {L6) 

where 

qn = qPh and q L = qP R . 

The validity of Eq. (1131) is general and does not refer to "massless" quark fields. 

From a phenomenological point of view the u and d quarks and to a lesser extent also the s quark 
have relatively small masses in comparison to a typical hadronic scale of the order of 1 GeV. On the other 
hand, we will neglect the three heavy quarks c, b, and t, because we will restrict ourselves to energies 
well below the production threshold of particles containing a heavy (anti-) quark. In the following, 
we will approximate the full QCD Lagrangian by its light-flavor version, and will consider the chiral 
limit for the three light quarks u, d, and s. To that end, we apply Eq. (FlBl) to the term containing the 
contraction of the covariant derivative with 7^. This quadratic quark form decouples into the sum of 
two terms which connect only left-handed with left-handed and right-handed with right-handed quark 
fields. The QCD Lagrangian in the chiral limit can then be written as 

£°QCD= E <1R,1 + QL,ti$ QL,l) - \$a^ U - (14) 

l=u,d,s 
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Note that because of Eq. ffTB"]) the quark-mass term generates a coupling between left- and right-handed 
quark fields. 



2.1.3 Global symmetry currents of the light quark sector 

Due to the flavor independence of the covariant derivative, £q CD is invariant under the infinitesimal 
global transformations of the left- and right-handed quark fields, 



Qr = 



(Z) - (l-ilXy-fc*)?*. (15) 

Note that the Gell-Mann matrices act in flavor space. £q CD is said to have a classical global U(3) L x 
U(3) R symmetry. Applying Noether's theorem Noether, 1918] , |Hill, 1951 , Gell-Mann and Levy, 1960 1 , 
from such an invariance one would expect a total of 2 x (8 + 1) = 18 conserved currents: 

LZ = q Ll ^q L , L» = q L rq L , K = fefy Qr, R" = QrI^Qr- (16) 



Making use of 

P l YPr ± PrYPl 
we introduce the linear combinations 



A£ = R£-I£ = q>f>y s ?±q, (18) 

which under a parity transformation of the quark fields, q(t,x) i— > 7og(i, —x), transform as vector and 
axial-vector current densities, respectively, 

P:V»(t,x)^V afM (t,-x), (19) 
P:A%(t,x)^-A afl (t,-x). (20) 

The conserved singlet vector current results from a transformation of all left-handed and right-handed 
quark fields by the same phase, 

V" = R» + = q^q. (21) 

The singlet axial-vector current originates from a transformation of all left-handed quark fields with 
one phase and all right-handed with the opposite phase, 

A» = - If = q^q. (22) 

Quantum fluctuations destroy the singlet axial- vector current conservation and there will be extra terms, 
referred to as anomalies |Bell and Jackiw, 1969 , Adler, 1969 , Adler and Bardeen, 1969 , resulting in 

dpA? = y^w^rer, e 0123 = i. 

The factor of three originates from the number of flavors. In the large N c (number of colors) limit of 
Ref. 't Hooft, 1974 the singlet axial-vector current is conserved, because the strong coupling constant 
behaves as g\ ~ iV" 1 . 
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2.1.4 Chiral algebra 

The invariance of £q CD under global SU(3) L x SU(3) fl x U(l) y transformations implies that also the 
QCD Hamilton operator in the chiral limit, Hq CB , exhibits a global SU(3) L x SU(3) H x U(l) y symmetry. 
As usual, the charge operators are defined as the space integrals of the charge densities, 

Q aL (t) = Jd 3 xq L (t,x)^-q L (t,x), (23) 

Qa R (t) = |d 3 2^(t,f)yrf,f), (24) 

Qvif) = Jd 3 x [q ] L (t } x)q L (t,x) + q R (t, x)q R (t, x)} . (25) 

For conserved symmetry currents, these operators are time independent, i.e., they commute with the 
Hamiltonian, 

[QaL, #QCd] = [QclR,Hq CD ] = [Q v ,Hq CD ] = 0. (26) 
The commutation relations among the charge operators reflect the underlying Lie algebra of SU(3) L x 
SU(3) fl x U(1) VI 

[QaL,QbL\ = ifabcQcL, (27) 
[QaR, QbR] = ifabcQcR, (28) 
[QaL,Qb R ] = 0, (29) 
[QaL,Q V ] = [QaR,Qv]=0. (30) 

Equations (1271) - (1301) are verified by expressing the commutators in terms of equal-time anti-commutation 
relations of the quark fields. 

It should be stressed that, even without being able to explicitly solve the equation of motion of the 
quark fields entering the charge operators of Eqs. (127]) - fl30l) . we know from the equal-time commuta- 
tion relations and the symmetry of the Lagrangian that these charge operators are the generators of 
infinitesimal transformations of the Hilbert space associated with Hq CD . Furthermore, their commu- 
tation relations with a given operator specify the transformation behavior of the operator in question 
under the group SU(3) L x SU(3) R x U(l) v . 

2.1.5 Quark masses and chiral symmetry breaking 

So far, we have discussed an idealized world with massless light quarks. The finite u-, d-, and s-quark 
masses explicitly break the chiral symmetry and generate divergences of the symmetry currents. As a 
consequence, the charge operators are, in general, no longer time independent. However, as first pointed 
out by Gell-Mann Gell-Mann, 1962] , the equal-time commutation relations still play an important role 



even if the symmetry is explicitly broken. 
Defining the quark-mass matrix as 

M = diag(m u , m d ,m s ), 

the quark-mass term in the QCD Lagrangian leads to a mixing of left- and right-handed fields [see 
Eq. (HI], 

C M = -qMq = -(q R Mq L + q L Mq R ). (31) 

Inserting the transformations of Eqs. (fT5j) into the quark-mass term of Eq. (13~T!) results in the variation 
SCm, from which one obtains for the divergences 



^ a = ~d^~ = ~ l 1 qL ~2 Mqn ~ QR M Y qL ) ' 

= —^ = -i(q L Mq R -q R Mq L ). (32) 
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The analogous expressions for d^R^ and <9 M i? M are obtained from Eqs. fl32l) through the substitution 
R <-> L. The divergences are proportional to the mass parameters which is the origin of the expression 
current-quark mass. In terms of the vector and axial-vector currents the divergences read 

d^ a = iq{M,^} l5 q, 
= 0, 

3o 2 

d»A» = 2iqMj 6 q + ^e ia/pa g^gSr, e 0123 = 1. (33) 

We are now in the position to summarize the various (approximate) symmetries of the strong inter- 
actions in combination with the corresponding currents and their divergences. 



In the limit of massless quarks, the sixteen currents and R£ or, alternatively, VJf and A% are 
conserved. The same is true for the singlet vector current V 1 , whereas the singlet axial- vector 
current A^ 1 has an anomaly. 

For any value of quark masses, the individual flavor currents u'j^u, d'j^d, and s7 M s are always 
conserved in the strong interactions reflecting the flavor independence of the strong coupling and 
the diagonality of the quark-mass matrix. Of course, the singlet vector current V^, being the sum 
of the three flavor currents, is always conserved. 

In addition to the anomaly, the singlet axial-vector current has an explicit divergence due to the 
quark masses. 

For equal quark masses, m u = = m s , the eight vector currents are conserved, because 
[A , 1] = 0. Such a scenario is the origin of the SU(3) symmetry originally proposed by Gell- 
Mann and Ne'eman Gell-Mann and Ne'eman, 1964| . The eight axial- vector currents A£ are not 



conserved. The divergences of the octet axial-vector currents of Eq. (|33|) are proportional to 
pseudoscalar quadratic forms. This can be interpreted as the microscopic origin of the PCAC 
relation (partially conserved axial- vector current) |Gell-Mann, 1964] , |Adler and Dashen, 1968 



which states that the divergences of the axial-vector currents are proportional to renormalized 
field operators representing the lowest-lying pseudoscalar octet. 

Taking m u = ma ^ m s reduces SU(3) flavor symmetry to SU(2) isospin symmetry. 

Taking m u ^ leads to isospin symmetry breaking. 



Various symmetry-breaking patterns are discussed in great detail in Ref. Pagels, 1975 



2.1.6 Green functions, chiral Ward identities, and generating functional 

For conserved currents, the spatial integrals of the charge densities are time independent, i.e., in a 
quantized theory the corresponding charge operators commute with the Hamilton operator. These 
operators are generators of infinitesimal transformations on the Hilbert space of the theory. The mass 
eigenstates should organize themselves in degenerate multiplets with dimensionalities corresponding to 
irreducible representations of the Lie group in question. For the moment, we assume that the dynamical 
system described by the Hamiltonian does not lead to a spontaneous symmetry breakdown. We will 
come back to this point later. Which irreducible representations ultimately appear, and what the actual 
energy eigenvalues are, is determined by the dynamics of the Hamiltonian. For example, SU(2) isospin 
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symmetry of the strong interactions reflects itself in degenerate SU(2) multiplets such as the nucleon 
doublet, the pion triplet, and so on. Ultimately, the actual masses of the nucleon and the pion should 
follow from QCD. 

It is also well-known that symmetries imply relations between S-matrix elements. For example, ap- 
plying the Wigner-Eckart theorem to pion-nucleon scattering, assuming the strong-interaction Hamilto- 
nian to be an isoscalar, it is sufficient to consider two isospin amplitudes describing transitions between 
states of total isospin / = 1/2 or I = 3/2. All the dynamical information is contained in these isospin 
amplitudes and the results for physical processes can be expressed in terms of these amplitudes together 
with geometrical coefficients, namely, the Clebsch-Gordan coefficients. 

In quantum field theory, the objects of interest are the Green functions which are vacuum expectation 
values of time-ordered products. Later on, we will also refer to matrix elements of time-ordered products 
between states other than the vacuum as Green functions. The physical scattering amplitudes are 
obtained from the Green functions using the reduction formalism |Lehmann et al., 1955| . Symmetries 
provide strong constraints not only for scattering amplitudes, i.e. their transformation behavior, but, 
more generally speaking, also for Green functions and, in particular, among Green functions. Even if 
a symmetry is broken, i.e., the infinitesimal generators are time dependent, conditions related to the 
symmetry-breaking terms can still be obtained using equal-time commutation relations. 

The symmetry currents relevant to the global SU(3) L x SU(3) R xU(l) y of QCD are given in Eqs. (TT7|) . 
( TTSj) . and (I2"TI) . Moreover, since we also want to discuss explicit symmetry breaking, we introduce the 
scalar and pseudoscalar densities 

S a = qX a q, P a = iqjsKq, a = 0, • • ■ , 8, (34) 

where Ao = y2/31. For example, linear combinations of S a and P a are needed to describe the diver- 
gences of the currents in Eqs. ( 1331) . Whenever it is more convenient, we will also use 

S(x) = q(x)q(x), P(x) = iq{x)^$q{x) , (35) 

instead of S and P . 

For example, the following Green functions of the "vacuum" sector 

(0\T[A^x)P b (y)}\0), 
(0\T[P a (x)J^y)P c (z)}\0), 
(0\T[P a (w)P b (x)P c (y)P d (z)]\0) 

are related to pion decay, the pion electromagnetic form factor (J^ 1 is the electromagnetic current), 
and pion-pion scattering, respectively. One may also consider similar time-ordered products evaluated 
between a single nucleon in the initial and final states in addition to the vacuum Green functions. This 
allows one to discuss properties of the nucleon as well as dynamical processes involving a single nucleon, 
such as 

(N\ J^(x)\N) <-> nucleon electromagnetic form factors, 

(N\A%(x)\N) <->• axial form factor + induced pseudoscalar form factor, 

(N\T[J^(x)J v (y)]\N) <-> Compton scattering, 

(iV|T[J M (x)P a (?/)]|iV) <-> pion electroproduction. 

Generally speaking, a chiral Ward identity relates the divergence of a Green function containing at 
least one factor of V£ or to some linear combination of other Green functions. The terminology chiral 
refers to the underlying SU(3) L x SU(3) R group. To make this statement more precise, let us consider 
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as a simple example the two-point Green function involving an axial-vector current and a pseudoscalar 
density, 



G' AP , ab (x,y) = (0\T[AZ(x)P b (y))\0) = B(x - y )(0\A»(x)P b (y)\Q) + e(y - x )(0\P b (y)A^x)\0), (36) 

and evaluate the divergence 

dlG^MiV) = d;[&(x - y )(0\AZ(x)P b (y)\0) + &(y - x )(0\P b (y)A»(x)\0)] 
= 5(x - y )(0\A a (x)P b (y)\0) - 5(x - y )(0\P b (y)A a (x)\0) 

+e(x - y )(0\d;A^x)P b (y)\0) + Q(y - x o ){0\P b (y)^ a (x)\0) 
= 5(xo-i/o)(0|[A«(x),P 6 ( 2/ )]|0) + (0|r[9XWn(l/)]|0), 

where we made use of d*Q(xo — yo) = 5(x — yo)gofj, = —d*Q(yo — x ). This simple example already 
shows the main features of (chiral) Ward identities. From the differentiation of the theta functions one 
obtains equal-time commutators between a charge density and the remaining quadratic forms. The 
results of such commutators are a reflection of the underlying symmetry. As a second term, one obtains 
the divergence of the current operator in question. If the symmetry is perfect, such terms vanish 
identically. If the symmetry is only approximate, an additional term involving the symmetry breaking 
appears. For a soft breaking such a divergence can be treated as a perturbation. 

The time ordering of n + 1 points x,xi, • ■ ■ , x n gives rise to (n+ 1)! distinct orderings, each involving 
products of n theta functions. Via induction, the generalization of the above simple example to an 
(n + l)-point Green function is symbolically of the form 

d*(0\T{r(x)A 1 (x 1 )---A n (x n )}\0) = 
(0\T{[d*J»(x)]A 1 (x 1 )---A n (x n )}\0) 
+5(x° - xl)(0\T{[J o (x), A 1 (x 1 )]A 2 (x 2 ) ■ ■ ■ A n (x n )}\0) 
+5(x° - x° 2 )(0\T{A 1 (x 1 )[J Q (x), A 2 (x 2 )} ■ ■ ■ A n (x n )}\0) 

+ ■■■ + 5(x° - x° n )(0\T{A 1 (x 1 ) ■ ■ ■ [J (x), A n (x n )}}\0), (37) 

where J M stands generically for any of the Noether currents. 

The discussion so far assumes that one explicitly works out the particular chiral Ward identity one 
is interested in. However, there is an elegant way of obtaining all chiral Ward identities from a single 
expression. To that end we introduce into the Lagrangian of QCD the couplings of the nine vector 
currents, eight axial-vector currents, nine scalar quark densities, and nine pseudoscalar quark densities 
to external c- number fields Gasser and Leutwyler, 19 85 1: 



£ — £qcd + ^cxt) (3£ 



where 



8 A 1 8 A 8 8 

a=l A 6 a=l A a=0 a=0 

= ^(^+3^) + l5^)q - q(s - i-f 5 p)q. (39) 

The 35 real functions v%(x), v? s -,(x), a^(x), s a (x), and p a {x), will collectively be denoted by [v, a, s,p). A 
precursor of this method was already used by Bell and Jackiw |Bell and Jackiw, 1969 in their discussion 



of the anomalous divergences in the ir° — > 77 decay. The Green functions of the vacuum sector may be 
combined in the generating functional 



exp(iZ[v , a, s,p\) = (0|Texp 



J d 4 xC cxt (x) 



|0) . (40) 
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Note that both the quark field operators q in £ ext and the ground state |0) refer to the chiral limit, 
indicated by the subscript in Eq. (T4U1) . A particular Green function is then obtained through a 
functional derivative with respect to the external fields. As an example, suppose we are interested in 
the scalar w-quark condensate in the chiral limit, (0|m , u|0)o. We express uu as 



uu 




and obtain 



(0|«(ar)tt(a;)|0)o 



5 



+ 



S 



+ 



s 



3Ss (x) Ss 3 (x) ^5s 8 (x) 



exp(iZ[v, a, s,p]) 



=p=0 



From the generating functional, we can even obtain Green functions of the "real world," where 
the quark fields and the ground state are those with finite quark masses. For example, the two-point 
function of two axial- vector currents of the "real world," i.e., for s = di&g(m u ,m d ,m s ), and the "true 
vacuum" |0), is given by 



(0|T[^( a? K(0)]|0>= (-if 



s 2 



6a a)1 (x)5a bl/ (0) 



exp(iZ[v, a, s,p]) 



ti=a=p=0,s=diag(m li ,,m^,m s ) 



(41) 



Note that the left-hand side involves the quark fields and the ground state of the "real world," whereas 
the right-hand side is the generating functional defined in terms of the quark fields and the ground state 
of the chiral limit. The actual value of the generating functional for a given configuration of external 
fields v, a, s, and p reflects the dynamics generated by the QCD Lagrangian. 

The (infinite) set of all chiral Ward identities resides in an invariance of the generating functional 
under a local transformation of the external fields Gasser and Leutwyler, 1984] , |Leutwyler, 1994 . The 
use of local transformations allows one to also consider divergences of Green functions. We require 
L of Eq. (1551) to be a Hermitian Lorentz scalar, to be even under P, C, and T, and to be invariant 
under local chiral transformations. In fact, it is sufficient to consider P and C, only, because T is then 
automatically incorporated owing to the CPT theorem. 

Under parity, the quark fields transform as 



p 

q f (t,x) i ► j q f (t,-x), 
and the requirement of parity conservation, 



C(t,x) £ C(t, 



-x) 



leads, using the results of Table HJ to the following constraints for the external fields, 



P 



■I'-- 



P (s) 



P 



cr h-> -a 



P 

s i— > s, 



p 

p I — > —p. 



(42) 



(43) 



(44) 



In Eq. (jHJ) it is understood that the arguments change from (t,x) to (t, —x). 
Similarly, under charge conjugation the quark fields transform as 



where the subscripts a and (3 are Dirac spinor indices, 



C 



(45) 



C 



z 7 y 



(° 








-1 


\ 








1 










-1 
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a*" 




7^75 


7or7o 
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7a. 




-7s 


-7a«75 



Table 4: Transformation properties of the Dirac matrices F under parity. 



r 


1 






75 


7 M 75 


—CT T C 


1 






75 


7 M 75 



Table 5: Transformation properties of the Dirac matrices T under charge conjugation. 

is the usual charge conjugation matrix, and / refers to flavor. Taking Fermi statistics into account, one 
obtains 

qTFq = -qCT T CF T q } 

where F denotes a matrix in flavor space. In combination with Table [5] it is straightforward to show 
that invariance of C ext under charge conjugation requires the transformation properties 

C T (s) G l s )T C T C T T t A a\ 

where the transposition refers to the flavor space. 

Finally, we need to discuss the requirements to be met by the external fields under local SU(3) L x 
SU(3) R x U(l)y transformations. In a first step, we write Eq. (139]) in terms of the left- and right-handed 
quark fields. Using the projection operators of Eq. ffTTj) the Lagrangian of Eq. (1391) reads 

£ = ^QCD + 9l7 M (l» + QL + q R ^ (?> + ^ s) ) q R - q R (s + ip)q L - q L (s - ip)q R . (47) 

Equation (j4"T|) remains invariant under local transformations 

q R h-> exp I -i— — I V R {x)q R , 

q L ^ exp (^-i^Y^J Vi{x)q L) (48) 

where V R (x) and Vl{x) are independent space-time-dependent SU(3) matrices, provided the external 
fields are subject to the transformations 



in >-> VLipVl-idpViyl 

s + ip h-> V R (s + ip)Vl, 

s-ip i-> Vl{s- ip)V ] R . (49) 

The derivative terms in Eq. (T4"9"j) serve the same purpose as in the construction of gauge theories, i.e., 
they cancel analogous terms originating from the kinetic part of the quark Lagrangian. Note that the 
external currents are coupled with an "opposite" sign in comparison with our convention for gauge 
theories. 



13 



There is another, yet, more practical aspect of the local invariance, namely: such a procedure allows 
one to also discuss a coupling to external gauge fields in the transition to the effective theory to be 
discussed later. For example, a coupling of the electromagnetic field to point-like fundamental particles 
results from gauging a U(l) symmetry. Here, the corresponding U(l) group is to be understood as a 
subgroup of a local SU(3) L x 811(3)^. Another example deals with the interaction of the light quarks 
with the charged and neutral gauge bosons of the weak interactions. 

Let us consider both examples explicitly. The coupling of quarks to an external electromagnetic 
field Afj, is given by 

r M = 1(1 = -eQA^, (50) 
where Q = diag(2/3, —1/3,-1/3) is the quark charge matrix and e > the elementary charge: 

/2_ 1 - 1 _ 

£ ext = -eA^{q L Q^q L + q R Q^q R ) = -eA^gQ'fq = -eA^ ( -uj»u - -dj^d - -S'fs 

On the other hand, if one considers only the two-flavor version of QCD one has to insert for the external 
fields 

ffi — l/j, — ~ e ~^~-At> v \? = ~-^A(t. (51) 

In the description of semi-leptonic interactions such as 7r~ — > (JTv^, tt~ — > ir°e~v e , or neutron 
decay n — > pe~u e one needs the interaction of quarks with the massive charged weak bosons = 

(Wi M TOV 2M )/V2, 

^ = 0, Z M =--^(V^T+ + H.c.) ) (52) 
where H.c. refers to the Hermitian conjugate and 

o v ud v us 
:/+ lo o o 



Here, Vij denote the elements of the Cabibbo-Kobayashi-Maskawa quark-mixing matrix describing the 
transformation between the mass eigenstates of QCD and the weak eigenstates Amsler et al., 2008 , 



\V ud \ = 0.97418 ± 0.00027, \V US \ = 0.2255 ± 0.0019. 

At lowest order in perturbation theory, the Fermi constant is related to the gauge coupling g and the 
W mass as 

a 2 

Gp = V^-^- = 1.16637(1) x 10" 5 GeV~ 2 . (53) 
8M W 



Making use of 

q L rW+T + q L = ^W+[V ud ur(l ~ 75)^ + V us uY(l - 7s)a], 

we see that inserting Eq. fl52|) into Eq. fT4T|) leads to the standard charged- current weak interaction in 
the light-quark sector, 

{W+iVutuYil - ls)d + V us uY(l - 75>] + H.c.} . 



--ext 



2V2 



The situation is slightly different for the neutral weak interaction. Here, the three-flavor version 
requires a coupling to the singlet axial- vector current which, because of the anomaly of Eq. (1331 . we 
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have dropped from our discussion. On the other hand, in the two-flavor version the axial-vector current 
part is traceless and we have 



,(«) 



,73 
2 

9 r 3 



etan(8 w )—Z fl , 



T3 



cos(^) 2 ^ + etan(My^, 
etan(6 l vi/) 



2, 



(54) 



where 9w is the weak angle. With these external fields, we obtain the standard weak neutral-current 
interaction 



9 



ext 



1 4 



2 cos^^vk, 
where we made use of e = g sin (6*^ 



- - - sm {V W 'j 



-- + -sin^(^; 



1 + <T 5 H' 



2.1.7 PCAC in the presence of an external electromagnetic field 

Finally, the technique of coupling the QCD Lagrangian to external fields also allows us to determine 
the current divergences for rigid external fields, i.e., fields which are not simultaneously transformed. 
For the sake of simplicity we restrict ourselves to the two-flavor sector. (The generalization to the 
three-flavor case is straightforward.) 

Consider a global chiral transformation only and assume that the external fields are not simultane- 



ously transformed. In this case the divergences of the currents read Fuchs and Scherer, 2003 



T<i T' T' T' 

«W[-o , vM + *ff7 M 75[-j, a^q - iq[-±, s]q - qj 5 [-^,p]q, 

7~" T' 7~" T' 

ni^lA^ V M + Wl*^ o M ] Q + nib{^,s}q + q{^,p}q. 



(55) 
(56) 
m u = m d 

in combination with a coupling to an external electromagnetic field [see Eq. ( loTT) . — = p]. 
The expressions for the divergence of the vector and axial-vector currents, respectively, are given by 



As an example, let us consider the QCD Lagrangian for a finite light quark mass rh 



Fuchs and Scherer, 2003 



-^eA^^q 

T, 



'^3ij^A^Vj , 



Ti 



eA^q^^—q + 2miq-f 5 —q = -eA^A^ + mP, 



(57) 
(58) 



with the isovector pseudoscalar density Pi = iqj^Tiq. In fact, Eq. ( 158]) is incomplete, because the third 
component of the axial- vector current, A%, has an anomaly which is related to the decay tt° — > 77. The 
full equation reads 



m 



P - eA^eujA^ + 6,, 



(3 



32tt 



2 ^-/lupa 



^0123 



(59) 



where = d^A v — d u A^ is the electromagnetic field strength tensor. 

We emphasize the formal similarity of Eq. fl58|) to the (pre-QCD) PCAC (Partially Conserved Axial- 
Vector Current) relation obtained by Adler |Adler, 1965] through the inclusion of the electromagnetic 
interactions with minimal electromagnetic coupling. Since in QCD the quarks are taken as truly ele- 
mentary, their interaction with an (external) electromagnetic field is of such a minimal type. In Adler's 
version, the right-hand side of Eq. (1591) contains a renormalized field operator creating and destroy- 
ing pions instead of rhPi. From a modern point of view, the combination mPij '(M%F n ) serves as an 
interpolating pion field. Furthermore, the anomaly term is not yet present in Ref. |Adler, 1965 . 
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2.2 Spontaneous symmetry breaking and Goldstone theorem 
2.2.1 Linear sigma model 

Spontaneous symmetry breaking occurs if the ground state has a lower symmetry than the Hamiltonian. 
For example, in the linear sigma model |Schwinger, 1957| , |Gell-Mann and Levy, 1960 , the Lagrangian 
is constructed in terms of the 0(4) multiplet (a, ir%, n 2 , 7r 3 ), 



£ = - (d.ad^a + ■ d^) -^(cr 2 + vr 2 



A 
4< 



7T 



=;2\2 



(60) 



where A > 0. Under parity we assume a(t,x) i— »■ a(t, —x) and iTi(t,x) h- > — n^t, —x). The Lagrangian is 
invariant under the infinitesimal transformations 
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The linear combinations Ri = (T, + Ti + %)/2 and Li = (Tj — Ti + %)/2, i = 1,2, 3, satisfy the commutation 
relations corresponding to an SU(2)xSU(2) Lie group. The multiplet (a, 7Ti, 7r 2 , 7r 3 ) transforms according 
to the (|, |) representation. By choosing m 2 < 0, the symmetry is realized in the Nambu- Goldstone 
mode |Nambu, i 960 



Goldstone, 1961 . Let us assume that the ground state is characterized by the 



vacuum expectation values 




(vr,) = 0. 



Introducing a 
1 



v + a' , the Lagrangian reads 
1 



C = - (<V<9V + <9 M 7? • d tM n) + -{-2m 2 )a' 2 + Xva'{7r 2 l + 7T 2 2 +7T 



>2\ i ^ /_2 i 2, 2, /2\2 ^4 
) + ^K+ 7r 2+ 7r 3+ (T ) • 



(61) 



(62) 



The ground-state configuration is no longer invariant under the full group 0(4). While the generators 
Tx, T 2 , and T 3 annihilate the ground state of Eq. (IBTi) . the generators T 4 , T 5 , and T 6 do not. The model- 
independent feature of the above example is given by the fact that for each of the three generators T 4 , T 5 , 
and Tq which do not annihilate the ground state one obtains a massless Goldstone boson. This is why 
Eq. (1621) contains no mass terms for the pions. In fact, the number of Goldstone bosons is determined 
by the structure of the symmetry groups |Goldstone et al., 1962| . Let G denote the symmetry group 
of the Lagrangian with nc generators and the subgroup H the symmetry group of the ground state 
with riH generators. For each generator which does not annihilate the vacuum one obtains a massless 
Goldstone boson, i.e., the total number of Goldstone bosons equals ric — nu- 

The Lagrangians used in motivating the phenomenon of a spontaneous symmetry breakdown are 
typically constructed in such a fashion that the degeneracy of the ground states is built into the potential 



16 



at the classical level (the prototype being the "Mexican hat" potential). As in the above case, it is then 
argued that an elementary Hermitian field of a multiplet transforming non-trivially under the symmetry 
group G acquires a vacuum expectation value signaling a spontaneous symmetry breakdown. However, 
there also exist theories such as QCD where one cannot infer from inspection of the Lagrangian whether 
the theory exhibits spontaneous symmetry breaking. Rather, the criterion for spontaneous symmetry 
breaking is a non-vanishing vacuum expectation value of some Hermitian operator, not an elementary 
field, which emerges through the dynamics of the underlying theory. In particular, we will see that the 
quantities developing a vacuum expectation value may also be local Hermitian operators composed of 
more fundamental degrees of freedom of the theory. 

While the model of Eq. fl60l) is constructed to illustrate the concept of a spontaneous symmetry 
breaking, it is not fully understood theoretically why QCD should exhibit this phenomenon. We will first 
motivate why experimental input, the hadron spectrum of the "real" world, indicates that spontaneous 
symmetry breaking happens in QCD. Secondly, we will show that a non-vanishing singlet scalar quark 
condensate is a sufficient condition for a spontaneous symmetry breaking in QCD. 

2.2.2 The hadron spectrum 

We saw in Section [2.1.31 that the QCD Lagrangian possesses an SU(3) L x SU(3) fl x U(l)y symmetry 
in the chiral limit in which the light quark masses vanish. From symmetry considerations involving the 
Hamiltonian Hq CB om y> one would naively expect that hadrons organize themselves into approximately 
degenerate multiplets fitting the dimensionalities of irreducible representations of the group SU(3) L x 
SU(3) H xU(l)y. The U(l) y symmetry results in baryon number conservation and leads to a classification 
of hadrons into mesons (B = 0) and baryons (B = 1). The linear combinations Q a y = Q a R + QaL and 
QaA = QaR. — QaL of the left- and right-handed charge operators commute with Hq CD , have opposite 
parity, and thus for states of positive parity one would expect the existence of degenerate states of 
negative parity (parity doubling) which can be seen as follows. 

Let \a, +) denote an eigenstate of Hq CB an d parity with eigenvalues E a and +1, respectively, 

#QCdI">+) = Ea\0t,+), 

P\a,+) = \a,+), 

such as, e.g., a member of the ground-state baryon octet (in the chiral limit). Defining \4> aa ) = 
Qo,a\®, +), because of [Hq CD , Q aA ] = 0, we have 

-^qcdI<?W = HQ CD Q aA \a, +) = Q a AHQ CD \a, +) = E a Q aA \a, +) = E a \(j) aa ), 
P\<p aa ) = PQaAP- 1 P\a,+) = -Q aA (+\a,+)) = -\<j )aa ). 

The state \4> aa ) can be expanded in terms of the members of a multiplet with negative parity, 

\<fiaa) = QaA\a,+) = \f3, -) (f3, - | Q aA | «, +) = t af3a \(5 } -). 

However, the low-energy spectrum of baryons does not contain a degenerate baryon octet of negative 
parity. Naturally the question arises whether the above chain of arguments is incomplete. Indeed, we 
have tacitly assumed that the ground state of QCD is annihilated by the generators Q aA . Let b^ a+ 
denote an operator creating quanta with the quantum numbers of the state Similarly, let 

create degenerate quanta of opposite parity. Expanding 

[QaA, = b^_t a f 3a , 

the usual chain of arguments then works as 

Qa A \a,+) = Q aA bl + \0) = (lQaA,bi + ] + bi + QoA)\0)=t a p a by\0). (63) 

^ 
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Figure 1: Pseudoscalar meson octet in an S) diagram. Baryon number B = 0. Masses in MeV. 



However, if the ground state is not annihilated by Q a A, the reasoning of Eq. fl63|) does no longer apply. In 
that case the ground state is not invariant under the full symmetry group of the Lagrangian resulting 
in a spontaneous symmetry breaking. In other words, the non-existence of degenerate multiplets of 
opposite parity points to the fact that SU(3) y instead of SU(3) £ x SU(3) fl is approximately realized as 
a symmetry of the hadrons. Furthermore, the octet of the pseudoscalar mesons is special in the sense 
that the masses of its members are small in comparison with the corresponding 1~ vector mesons. They 
are the candidates for the Goldstone bosons of a spontaneous symmetry breaking. 

According to the Coleman theorem Coleman, 1966 , the symmetry of the ground state determines 
the symmetry of the spectrum, i.e. 

Q aV \0) = Qv\0) = (64) 

implies SU(3) y multiplets which can be classified according to their baryon number. In the reverse 
conclusion, the symmetry of the ground state can be inferred from the symmetry of the spectrum. 
Figures CD and [2] show the octets of the lowest-lying pseudoscalar-meson states and the lowest-lying 
baryon states of spin-parity | + , respectively. 

The axial charges satisfy the commutation relations 



[QaA, QbA\ 
[QaVi QbA] 



ifabcQcVi 
ifabcQcA- 



Since the parity doubling is not observed for the low-lying states, 
annihilate the ground state, 

QaA\0) ^ 0, 



(65) 
(66) 

one assumes that the Q a A do not 

(67) 



i.e., the ground state of QCD is not invariant under "axial" transformations. In the present case, 
G = SU(3) L x SU(3) R with n G = 16 and H = SU(3) y with n H = 8 and we expect eight Goldstone 
bosons. According to the Goldstone theorem |Goldstone, 1961] , |Goldstone et al., 1962| , to each axial 
generator Q a A, which does not annihilate the ground state, corresponds a massless Goldstone boson 
field (j) a (x) with spin 0, whose symmetry properties are tightly connected to the generator in question. 
The Goldstone bosons have the same transformation behavior under parity, 



<j>a{t,X) H-> ~4> a (t, -X), 



(6£ 
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Figure 2: Baryon octet (J = |) in an (J3, 5) diagram. Masses in MeV. Baryon number B = 1. 

i.e., they are pseudoscalars, and transform under the subgroup H = SU(3) y , which leaves the vacuum 
invariant, as an octet [see Eq. (|66l)]: 



[QaV,(pb(x)} = ifabc<Pc{x). 



(69) 



2.2.3 The scalar singlet quark condensate 



In the following, we will show that a non-vanishing scalar singlet quark condensate in the chiral limit 
is a sufficient (but not a necessary) condition for a spontaneous symmetry breaking in QCD. In this 
section all physical quantities such as the ground state, the quark operators etc. are considered in the 
chiral limit. 

Let us first recall the definition of the nine scalar and pseudoscalar quark densities: 



S a (y) = q(y)Kq(y), a = 
P a (y) = iq(y)i5Kq(y), a = o 



1 1 w 5 



(70) 
(71) 



where Ao = y 2/3 1. We need the equal-time commutation relation of two quark operators of the form 

Ai(x) = q*(x)Aiq(x), where A4 symbolically denotes Dirac and flavor matrices and a summation over 
color indices is implied: 

[A^t, x), A 2 (t, y)} = 5 3 {x - y)q\x)[A u A 2 ]q{x). 



(72) 



With the definition 



and using 



Q aV (t) = I dW{t,xY-^q(t,x) 



[y, 7oA ] = 0, [ y , 7 A b ] = ^ifabcK 



we see, after integration of Eq. ( 1721) over x, that the scalar quark densities of Eq. ( 1701 transform under 
SU(3) V as a singlet and as an octet, respectively, 



[Qav{t),S (y)\ 



0, a 



(73) 
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[Qav{t), Sb 



Y.fabcS c (y), a,b= 8, 



(74) 



c=l 



with analogous results for the pseudoscalar quark densities. Using the relation 



fabcfabd — ^cd 

a, 6=1 



(75) 



for the structure constants of SU(3), we re-express the octet components of the scalar quark densities 



as 



S a (y) 



^ fabc[Qbv(t), S t 



(76) 



b,c=l 



In the chiral limit the ground state is necessarily invariant under SU(3)y Vafa and Witten, 1984 
i.e., Qavlfy — 0, and we obtain from Eq. (ITS"]) 



(0\S a (y)\0) = (0|5 a (0)|0) = (S a ) = 0, a = 1 



5 5 ^5 



(77) 



where we made use of translational invariance of the ground state. As an intermediate result we see 
that the octet components of the scalar quark condensate must vanish in the chiral limit. From Eq. 
f }77j) . we obtain for a = 3 

(uu) — (dd) = 0, 



i.e. (uu) = (dd) and for a = 8 



{uu) + (dd) - 2(ss) = 0, 



i.e. (uu) = (dd) = (ss). 

Because of Eq. (!73|) a similar argument cannot be used for the singlet condensate, and if we assume 
a non- vanishing singlet scalar quark condensate in the chiral limit, we find using (uu) = (dd) = (ss): 



7^ (<?<?) = ( uu + dd+ ss) = 3(uu) = 3(dd) = 3{ss). 
Finally, we make use of (no summation implied!) 



(78) 



W 2 [75y,7o75AJ = ^llo 



in combination with 



Xl — Xl — X 



to obtain 



/ 1 
1 

V o o o 



— \ 2 



/ 1 


V o o i 



Ag — Xj 




\ 2 



1 

3 



1 
1 
4 



i[QaA(t),P a (y)} 



(79) 



uu + dd, a = 1, 2, 3 

uu + ss, a = 4, 5 

dd + ss, a = 6, 7 

\(uu + dd + Ass), a = S 

where we have suppressed the y dependence on the right-hand side. We evaluate Eq. f!79|) for a ground 
state which is invariant under SU(3)y, assuming a non-vanishing singlet scalar quark condensate, 



\i[QaA(t),P a (y)]\0) = -(qq), a = 1, 



(80) 
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where, because of translational invariance, the right-hand side is independent of y. Inserting a complete 
set of states into the commutator of Eq. fISUl) yields that both the pseudoscalar density P a (y) as well 
as the axial charge operators Q a A must have a non-vanishing matrix element between the vacuum and 
massless one-particle states In particular, because of Lorentz covariance, the matrix element of the 
axial- vector current operator between the vacuum and these massless states, appropriately normalized, 
can be written as 

(OK(O)|0 6 (p)) = i^F o <U, (81) 

where Fq w 93 MeV denotes the "decay" constant of the Goldstone bosons in the chiral limit. From 
Eq. flHTl) we see that a non-zero value of Fq is a necessary and sufficient criterion for spontaneous chiral 
symmetry breaking. On the other hand, because of Eq. ( |80l) a non-vanishing scalar quark condensate 
(qq) is a sufficient (but not a necessary) condition for a spontaneous symmetry breakdown in QCD. 



3 Mesonic chiral perturbation theory 

Our goal is the construction of the most general theory describing the dynamics of the Goldstone 
bosons associated with the spontaneous symmetry breakdown in QCD. In the chiral limit, we want the 
effective Lagrangian to be invariant under G = SU(3) L x 811(3)^ x U(l) y . It should contain exactly eight 
pseudoscalar degrees of freedom transforming as an octet under the subgroup H = SXJ(3) V . Moreover, 
taking account of spontaneous symmetry breaking, the ground state should only be invariant under 
SU(3) V x U(l) v . 



3.1 Transformation properties of the Goldstone bosons 

The purpose of this section is to discuss the transformation properties of the field variables describing 



the Goldstone bosons Weinberg, 1968 , Coleman et al., 1969 , Callan et al., 19 69|. We will need the 



concept of a nonlinear realization of a group in addition to a representation of a group which one usually 
encounters in Physics. We will first discuss a few general group-theoretical properties before specializing 
to QCD. 



3.1.1 General considerations 

Let us consider a physical system described by a Lagrangian which is invariant under a compact Lie 
group G. We assume the ground state of the system to be invariant under only a subgroup H of G, 
giving rise to n = no — % Goldstone bosons. Each of these Goldstone bosons will be described by an 
independent field & which is a smooth real function on Minkowski space M 4 . These fields are collected 
in an n-component vector $ = (0i, • • • , (f) n ), defining the vector space M%. Our aim is to find a mapping 
tp which uniquely associates with each pair (g, $)gGx Mi an element ip(g, $) G Mi with the following 
properties: 

<p(e, $) = $V$€M lf e identity of G, (82) 
<p(gi,<p{g2,$)) = vigito,®) V 01,02 e G, V$ e Ml (83) 

Such a mapping defines an operation of the group G on M\. The construction proceeds as follows 
|Leutwyler, 19921 . Let $ = denote the ' 'origin" of M\ which, in a theory containing Goldstone bosons 
only, loosely speaking corresponds to the ground state configuration. Since the ground state is supposed 
to be invariant under the subgroup H we require the mapping p to be such that all elements h € H 
map the origin onto itself. In this context the subgroup H is also known as the little group of $ = 0. 

We will establish a connection between the Goldstone boson fields and the set of all left cosets 
{gH\g G G} which is also referred to as the quotient G/H. For a subgroup H of G the set gH = 
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{gh\h E H} defines the left coset of g (with an analogous definition for the right coset) which is one 
element of G/H. For our purposes we need the property that cosets either completely overlap or are 
completely disjoint, i.e, the quotient is a set whose elements themselves are sets of group elements, and 
these sets are completely disjoint. Under all elements of a given coset gH the origin is mapped onto 
the same vector in R n : 

<p(gh, 0) = (p(g, (p{h, 0)) = <p(g, 0) V g E G and h E H. 

Secondly, the mapping is injective with respect to the elements of G/H. Consider two elements g and 
g' of G where g' £ gH. Let us assume <p(g, 0) = (pig', 0): 

= (pie, 0) = <p(g- l g, 0) = (pig' 1 , (pig, 0)) = <p(g-\ (pig', 0)) = (pig^g', 0). 

However, this implies g~ 1 g' E H or g' E gH in contradiction to the assumption g 1 E" gH and therefore 
(pig, 0) = (pig 1 ', 0) cannot be true. In other words, the mapping can be inverted on the image of (pig, 0). 
The conclusion is that there exists an isomorphic mapping between the quotient G/H and and the 
Goldstone boson fields. Of course, the Goldstone boson fields are not constant vectors in R™ but 
functions on Minkowski space. This is accomplished by allowing the cosets gH to also depend on x. 

Now let us discuss the transformation behavior of the Goldstone boson fields under an arbitrary g E 
G in terms of the isomorphism established above. To each $ corresponds a coset gH with appropriate 
g. Let / = gh E gH denote a representative of this coset such that 

Q> = <p(f,0) = cp(gh,0). 

Now apply the mapping (pig) to <3>: 

(p{g, $) = (p{g, (p^h, 0)) = (pig~gh, 0) = <p{f, 0) = f G gi$H). 

In order to obtain the transformed $' from a given $ we simply need to multiply the left coset gH 
representing $ by g in order to obtain the new left coset representing 

I T 
~gH ± g~gH 

This procedure uniquely determines the transformation behavior of the Goldstone bosons up to an 
appropriate choice of variables parameterizing the elements of the quotient G/H. 



3.1.2 Application to QCD 

The symmetry groups relevant to the application in QCD are 

G = SU(iV) x SU(iV) = {(L, R)\L E SXJ(N),Re SXJ(N)} and H = {(V,V)\V E SU(iV)} = SU(iV). 
Let g = (L, R) E G. We characterize the left coset gH through the SU(iV) matrix U = RL* 



Bal achandran et al., 1991 such that gH — (1, RU)H. This corresponds to the convention of choosing 



as the representative of the coset the element which has the unit matrix in its first argument. The 
transformation behavior of U under g = (L, R) E G is obtained by multiplication in the left coset: 

ggH = (L, RB.D)H = (1, RRDl^)(L, L)H = (1, R{RL ] )L ] )H, 



i.e. 



U = RD \-*U' = R{RV)L ] = WIS. (84) 
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As mentioned above, we need to introduce an x dependence to account for the fact that we are dealing 
with fields: 

U(x) ^ RU(x)lS. (85) 
For the physically relevant cases the corresponding unitary matrices may be parameterized as 

U(x) = exp ( j , (86) 



where, for N = 2, 



and, for iV = 3, 



^ I (87) 



i=i 



v/2' 



7T — 7T 



8 
a=l 



/ 7T° + ^77 V / 27T+ \ 

V J 



The origin </>(x) = 0, i.e. [/q = 1, denotes the ground state of the system. Under transformations of 
the subgroup H = {{V,V)\V e SU(iV)} corresponding to rotating both left- and right-handed quark 
fields in QCD by the same V, the ground state remains invariant, 

U Q ^ VUqV^ = VV ] = 1 = U . 

On the other hand, under "axial transformations," i.e. rotating the left-handed quarks by A and the 
right-handed quarks by A\ the ground state does not remain invariant, 

U ^ A ] U Q A ] = A ] A ] ^ U , 

which is consistent with the assumed spontaneous symmetry breakdown. 

The traceless and Hermitian matrices of Eqs. (|87l) and (l8~8~I) contain the Goldstone boson fields. 
We want to discuss their transformation behavior under the subgroup H = {(V,V)\V G SU(iV)}. 
Expanding 

we immediately see that the transformation behavior of Eq. (j85p restricted to the subgroup H, 



F 2F 2 • ' v ~ ' ~F 2F 2 • ' ~ ' 2if 



1 + l ~ ~ ^2 + • ■ ■ ^ V(l + l— - — ^ + • • 0^ = 1 + + 



implies 

<f> ^ V(f)V ] . (89) 

However, this corresponds exactly to the fact that the Goldstone bosons transform according to the 
adjoint representation under SU(3)y, i.e. they transform as an octet. 

For group elements of G of the form (A, A^) one may proceed in a completely analogous fashion. 
However, one finds that the fields <p a do not have a simple transformation behavior under these group 
elements. 
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3.2 Effective Lagrangian and power-counting scheme 

The application of effective field theory (EFT) to strong interaction processes has become one of the most 
important theoretical tools in the low-energy regime. The basic idea consists of writing down the most 
general possible Lagrangian, including all terms consistent with assumed symmetry principles, and then 
calculating matrix elements with this Lagrangian within some perturbative scheme [Weinberg, 19 79 1. 
A successful application of this program thus requires two main ingredients: 

(1) A knowledge of the most general effective Lagrangian; 

(2) an expansion scheme for observables in terms of a consistent power-counting method. 
3.2.1 The lowest-order effective Lagrangian 

In terms of the SU(3) matrix U(x) of Eqs. 0861) and (I88p the most general, chirally invariant, effective 
Lagrangian with the minimal number of derivatives reads 

£ eff = ^Tr (dpUPrf) , (90) 

where Fq « 93 MeV is a free parameter which later on will be related to the pion decay 7r + — > 
H^Vfj.. Because of the trace property Ti(AB) = Tt(BA), the Lagrangian is invariant under the global 
SU(3) L x SU(3) fl transformation U \— > RJJU of Eq. (1851) . The global U(l)y invariance is trivially 
satisfied, because the Goldstone bosons have baryon number zero, thus transforming as i— > under 
U(l)y which also implies U i— > U. 

The substitution <f) a (t,x) i— > —<p a (t,x) or, equivalently, U(t,x) \— > W(t,x) provides a simple method 
of testing, whether an expression is of so-called even or odd intrinsic parity, i.e., even or odd in the 
number of Goldstone boson fields. For example, the Lagrangian of Eq. fl90|) is even. Since the Goldstone 
bosons of QCD are pseudoscalars, the true parity transformation is given by <fi a (t,x) i— > —<fi a (t, —x) or, 
equivalently, U(t,x) i— > W(t,—x). 

The purpose of the multiplicative constant Fq/A in Eq. (|9"U|) is to generate the standard form of 
the kinetic term ^d^ad^cpa, which can be seen by expanding the exponential U = 1 + i(fi/F + • • •, 
dfJJ = id^(p/F + • • •, resulting in 

+ ■ • • = ^d^ a d^ b Tr(X a X b ) + ■■■= X -d^ a d^ a + Ant, 

where we made use of Tr(A a Af,) = 2<5 afe . Since there are no other terms containing only two fields, the 
eight fields a indeed describe eight independent massless particles. 

There are no other independent, chirally invariant terms containing only two derivatives. A term of 
the type Tr[(<9 ^8^11)11^] may be re-expressed as 

Ti[(d^U)U ] ] = d^TriVUrf)) - Tr^l/dptf), 

i.e., up to a total derivative it is proportional to the Lagrangian of Eq. (!90j) . However, in the present 
context, total derivatives do not have a dynamical significance, i.e. they leave the equations of motion 
unchanged and can thus be dropped. The product of two invariant traces is excluded at lowest order, 
because Trfp^UW) = 0. 

Let us turn to the vector and axial-vector currents associated with the global SU(3) L x SU(3) H 
symmetry of the effective Lagrangian of Eq. (1901) . To that end, we parameterize 

L = l-^ L y, (91) 
R = t-ie^. (92) 
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In order to construct J% L , set = and choose = e^(x). Using 

U*U =1 d»{lPU) =0 d»U ] U = -U j d^U, (93) 

the variation of the Lagrangian can be brought into the form 

5C cS = ^f-id^Tr (X a d^u) . 



Applying the method of Ref. Gell-Mann and Levy, 1960 , we obtain for the left currents 



and, completely analogously, choosing e„ = and = (x), for the right currents 

j - = = -4 Tr ( wf/t ) • (95) 

Combining Eqs. ( 1941) and ( 195]) the vector and axial-vector currents read 

J'av = Jtt R + JZ L = -i^Tr(\ a [U,d^)), (96) 
Ka = J&-JL = -i^Tr(Xa{U,Zrrt}). (97) 

Furthermore, because of the symmetry of C e s under SU(3) L x SU^)^, both vector and axial-vector 
currents are conserved. Using the substitution U «-> ift and Eq. (1931) . the vector current densities J^ v 
of Eq. (|96|) contain only terms with an even number of Goldstone bosons and the axial-vector current 
densities J^ A of Eq. ( 19TI) only terms with an odd number of Goldstone bosons. To find the leading term 
let us expand Eq. ( 1971) in the fields, 

^=-4 Tr ( A 4 i+ -'-'^^ + --]) = ~ Fo ^ + ''-- 

We conclude that the axial- vector current has a non-vanishing matrix element when evaluated between 
the vacuum and a one-Goldstone boson state: 

(0| J^{x)\(j> b {p)) = (0| - Fod^M^lMp)) = -Fod^eM-iP ' = ^ M F expHp • x)5 ab . (98) 
Equation (I98p is the manifestation of Eq. (I8"T|) at lowest order in the effective field theory. 

3.2.2 Symmetry breaking through the quark masses 

Up to now, we have assumed a perfect SU(3) L x SU(3) i? symmetry. As has been discussed in Section 
12.1.51 the quark-mass term of QCD results in an explicit symmetry breaking, 

/ m u \ 

C M = -q R Mq L - q L M^q R , M=\ m d . (99) 

\ m s j 

In order to incorporate the consequences of Eq. (1991) into the effective-Lagrangian framework, one makes 
use of the following argument |Georgi, 1984| : Although M. is in reality just a constant matrix and does 
not transform along with the quark fields, Cm 01 Eq. (1991) would be invariant if Ai transformed as 



M ^ RMtf. (100) 
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One then constructs the most general Lagrangian C(U, M) which is invariant under Eqs. fl85l) and (11001) 
and expands this function in powers of A4. At lowest order in Ai one obtains 

£ s .b. = S£±Tr(Mrf + UM*), (101) 

where the subscript s.b. refers to symmetry breaking. In order to interpret the new parameter Bq, let 
us consider the Hamiltonian density corresponding to the sum of the Lagrangians of Eq. ( 190]) and (jTUTJ) : 

H cS = ^fTr(Utf) + ^Tr(W ■ W f ) - £ s . b . 

Since the first two terms are always larger or equal to zero, Ti, e e is minimized by constant and uniform 
fields. Using the ansatz 

= 00 + 02 + T^04 H 

for the minimizing field values and organizing the individual terms in powers of 1/F 2 , one finds = 
as the classical solution even in the presence of quark-mass terms. Now consider the energy density of 
the ground state (U min = U = t), 

(^cn)min = -F 2 B (m u + m d + m a ), (102) 

and compare its derivative with respect to (any of) the light quark masses m q with the corresponding 
quantity in QCD, 



«9(0|H QC d|0> 



= l(o\qq\o)o = l(m)o, 

m u =m c i=i r n s =0 



where (qq)o is the chiral quark condensate of Eq. f!78|) . Within the framework of the lowest-order 
effective Lagrangian, the constant Bq is thus related to the chiral quark condensate as 

3F Q 2 5 = -(qq} . (103) 

Let us add a few remarks. 

1. A term Ti(Ai) by itself is not invariant. 

2. The combination Tr(.Mt/^ — UA4<) has the wrong behavior under parity 0(t, x) — > — 0(t, — x). 

3. Because M. = , £ s .b. contains only terms even in 0. 

In order to determine the masses of the Goldstone bosons, we identify the terms of second order in the 
fields in £ s .b., 

£ s . h = -^Tr(^M) + ---. (104) 

For the sake of simplicity we consider the isospin-symmetric limit m u = rrid = rh so that the tt i] term 
vanishes and there is no TT°-r] mixing. The masses of the Goldstone bosons, to lowest order in the quark 
masses, are then given by 

Ml = 250m, (105) 
M 2 K = B (m + m s ), (106) 

Ml = ^B (m + 2m s ). (107) 
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Heisenberg ferromagnet 


QCD 


Symmetry of Hamiltonian 


0(3) 


SU(3) L x SU(3) i? 


Symmetry of 0) 


0(2) 


SU(3) y 


Vacuum expectation value 


(M> 


(qq)o 


Explicit symmetry breaking 


External magnetic field 


Quark masses 


Interaction 


— (A?) ■ H 


(H cS ) of Eq. (UU2J) 



Table 6: Comparison between the symmetry-breaking patterns of a Heisenberg ferromagnet and QCD. 



These results, in combination with Eq. (I103p . B = — (qq) Q / (3Fq), correspond to relations obtained 
in Ref. Gell-Mann et al., 1968] and are referred to as the Gell-Mann, Oakes, and Renner relations. 



Without additional input regarding the numerical value of B Q , Eqs. (I105P - (11071) do not allow for 
an extraction of the absolute values of the quark masses in and m s , because re-scaling Bq — > XBq in 
combination with m q — > m q /X leaves the relations invariant. For the ratio of the quark masses one 
obtains, using the empirical values M v = 135 MeV, = 496 MeV, and M v = 547 MeV, 

M 2 K _m + m s ^ !^_25 9 



Ml 2m m 



Wis 

—^ = 24.3. 108 



M% 2m s + m 
M% 3rh m 

Let us conclude this section with a remark on (qq)o- A non-vanishing quark condensate in the chiral 
limit is a sufficient but not a necessary condition for a spontaneous chiral symmetry breaking. The 
effective Lagrangian term of Eq. (HOip not only results in a shift of the vacuum energy but also in finite 
Goldstone boson masses and both effects are proportional to the parameter B . We recall that it was a 
symmetry argument which excluded a term Tr(A^) which, at leading order in A4, would decouple the 
vacuum energy shift from the Goldstone boson masses. The scenario underlying C s b. of Eq. (1 10 II) is 
similar to that of a Heisenberg ferromagnet which exhibits a spontaneous magnetization (M), breaking 
the 0(3) symmetry of the Heisenberg Hamiltonian down to 0(2). In the present case the analogue of 
the order parameter (M) is the quark condensate (qq)o- In the case of the ferromagnet, the interaction 
with an external magnetic field H is given by — (M) ■ H, which corresponds to Eq. (I102p . with the 
quark masses playing the role of the external field H (see Table [6]). However, in principle, it is also 
possible that B vanishes or is rather small. In such a case the quadratic masses of the Goldstone 
bosons might be dominated by terms which are nonlinear in the quark masses, i.e., by higher-order 
terms in the expansion of C(U,Ai). Such a scenario is the origin of the so-called generalized chiral 
perturbation theory |Knecht et al., 1995] . The analogue would be an antiferromagnet which shows a 
spontaneous symmetry breaking but with (M) = 0. The analysis of the s-wave 7T7r-scattering lengths 
|Colangelo et al., 2000] , Colangelo et al., 2001a] supports the conjecture that the quark condensate is 



indeed the leading order parameter of the spontaneously broken chiral symmetry (see also Sec. 13.2.41) . 
3.2.3 Construction of invariants 

So far, we have discussed the lowest-order effective Lagrangian for a global SU(3) L x SU(3) R symmetry. 
In Sec. I2.1.6l we stated that the Ward identities of QCD are obtained from a locally invariant generating 
functional involving a coupling to external fields. Therefore, following Refs. |Gasser and Leutwyler, 1984] , 



27 



|Gasser and Leutwyler, 1985 , we will promote the global symmetry of the effective Lagrangian to a local 
one, 

L — > V L (x), R — > V R (x), 

and introduce a coupling to the same external fields v, a, s, and p as in QCD [see Eq. (1471) ]. The 
transformation behavior of the special unitary matrix U of Eqs. (I86I) and (1581 under G = SU(3) L x 
SU(3) R , parity P, and charge conjugation C is given by 

U £ V R UVl, U(x, t) £ U\-x, t), U & U T . 

Given an object A transforming as VrAVI, such as U or Xi the covariant derivative of A, D^A, is 
defined as 

Dp A = d„A - ir^A + iAl„ h-> Z^A' = VrP^vJ. (109) 

The defining property is that the covariant derivative should transform as the object it acts on. In 
particular, the covariant derivative of U is given by 



D fl U = d fl U-tr fl U + iUl fl . 
For the external fields we introduce corresponding field strength tensors in matrix form as 



fa = <V<, - d » r » - r A & VnfSyi 
fa = dJv-dJ*-i[iM&VLfLvl. 



rBT/t 



(110) 

(111) 

(112) 



They are traceless, because Tr(/ M ) = Tr(r M ) = and the trace of any commutator vanishes. Fi- 
nally, we introduce the linear combination \ = 27?o(s + ip), where, e.g., pure QCD is given by 
X = 2B dmg(m u ,m d ,m s ). 

The effective Lagrangian is constructed in terms of U, W , x, x\ fuvi fuv an d covariant derivatives 
of these objects. Suppose we have matrices A, B, C, ■ ■ •, all of which transform as 

A&VrAVL b&v r bvI .... 

Invariants may be formed by "multiplying" in the following way: 

Tr(AB f ) £ Tr{V R A V ] L V L B ] V ] R ) = Tr(Vj V R AB ] ) = Tr(AB r ), 

1 

where the generalization to a longer string of terms is obvious and the product of invariant traces is 
also invariant: 

Tv(AB j CD^), Tr(AB^)Tr(CD^), (113) 
In the chiral counting scheme the elements count as 

17 = 0(5°), D fl U = 0(q), r^l^ = 0{q), fa R = 0(q 2 ), X = 0(q 2 ). 

Any additional covariant derivative counts as 0(q). The list of objects A up to and including order q 2 
which transform as A' = VrAVI reads 

U, D,U, D,D V U, X, Ufa, faU 

The construction of chirally invariant expressions up to and including order q 2 proceeds as follows. At 
O(q ) the only invariant term is a constant, Tr (yW) = Tr(l) = const. Because of Tr (d^UW) = 0, 
terms of the type Tr[C(g)] x Tr(- • •) are excluded. At 0(q 2 ) we have 



Tr (D^DvUrf) = -Tr [d v U{D^] , Tr [d^U^U^] , Tr [u(D^D u U) 
Tt ( X rt), Tr ([//), Tr[{Ufa)rf]=Tr(fa)=0, Tr (fa) = 0. 
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-Tr 



D^U{D U U) ] 



Because of Lorentz invariance, indices have to be contracted and the remaining three candidates are 



Tr 

Finally, due to parity conservation, 



D^UiD^Uy] , Tr (xU j ± Ux 



C(x, t) i — ► C{— x, t). 

Tr(x?7' — Ux ) has to be excluded because of the wrong parity. At 0(q 2 ), charge conjugation does not 
generate any additional constraint. 

The locally invariant lowest-order Lagrangian £ 2 is given by 



C 2 = ^Tr 



F 2 „ r i F? 



DpUiptUy] + -^Tr (xU ] + U X ] ) ■ (114) 

At 0(q 2 ) it contains two parameters: the SU(3) chiral limit of the Goldstone boson decay constant 
F » 93MeV, and, hidden in the definition of x, B = -(0|gg|0) /(3F 2 ). 

The lowest-order equation of motion corresponding to Eq. (11X41) is obtained by considering small 
variations of the SU(3) matrix, 



U'{x) = U(x) + 5U{x) = (l+*£ A a(^)Aal U(x), 



(115) 



where the A a (x) are real functions. The matrix U 1 satisfies both conditions U'U'* = 1 and det(U') = 1 
up to and including terms linear in A . Applying the principle of stationary action, the variation of the 
action reads 



SS = j 2 dt [ d 3 x ^ A a (x)Tr {\ a [D^UU ] - U{D^U) ] - X U ] + U X ] ]} , 
4 Jti J a=1 

where we made use of partial integration, the standard boundary conditions A a (t 1; x) = A a (t 2 , x) = 0, 
the divergence theorem, and the definition of the covariant derivative of Eq. f)109p . Since the test 
functions A a (x) may be chosen arbitrarily, we obtain eight Euler-Lagrange equations 

Tr{\ a [D 2 UlP-U(D 2 Uy -xrf + Urf]} = 0, a = l,---,8, (116) 
which may be combined into a compact matrix form 

( kL(U) = D 2 UU^ - U(D 2 UY - X U ] + U X ] + ^Tr( X C/t - U X r ) = 0. (117) 

The trace term in Eq. (11171) appears, because Eq. (11161) contains eight and not nine independent 
equations. 



3.2.4 Two simple applications: Pion decay and txtx scattering 

The Lagrangian £ 2 of Eq. (I114p has predictive power, once the low-energy coupling constant F is 
identified. This LEC may be obtained from the weak decay of the pion, 7r + — > fi + u^. For that purpose 
we insert the corresponding external fields of Eq. (|52|) . describing the interaction of quarks with the 
massive charged weak bosons, into £ 2 - The coupling of a single W boson to a single Goldstone boson 
originates from the covariant derivatives in £ 2 > 

^Tr[D,U(D»Uy] = i^Tril^tfU) + ■■■ = ^Tr(/ M ^0) + ■ • • , 
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and is given by 



C W4> = --^Tr[(W M + T + +W;T_)d>] = - g ^[W+(V ud d^+V us d^K-)+W;(V ud d^ + +V us d^K + )}. 

(118) 

The invariant amplitude of the weak pion decay is of the structure "leptonic vertex x W propagator x 
hadronic vertex," 



M 



9 



2V2 



]9pa_ 



-g^v ud (-ip 



-G F V ud F u u J(l - 75 K 



(119) 



where Gp is the Fermi constant of Eq. fl53|) and p denotes the four-momentum of the pion. In the 
gauge-boson propagator, momenta p have been neglected in comparison to the gauge-boson mass M w . 
The corresponding decay rate is 



r 2 T/ 2 

(j F v ud 172 2 



An 



The constant F Q is referred to as the pion-decay constant in the chiral limit. It measures the strength of 
the matrix element of the axial-vector current operator between a one-Goldstone-boson state and the 
vacuum [see Eq. (18~T1) ]. Since the interaction of the W boson with the quarks is of the V — A type and 
the vector current operator does not contribute to the matrix element between a single pion and the 
vacuum, pion decay is completely determined by the axial- vector current. The degeneracy of a single 
coupling constant F is removed at next-to-leading order, 0(q A ) |Gasser and Leutwyler, 1985 



once 



SU(3) symmetry breaking is taken into account. The empirical numbers for F w and Fk are 92.4 MeV 
and 113 MeV, respectively |Amsler et al., 
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Now that the LEC Fq has been identified, we will show how the lowest-order Lagrangian predicts the 
prototype of a Goldstone-boson reaction, namely, im scattering. We consider C 2 in the SU(2) L x SU(2) R 



sector with r„ = /, 



0. 



C 2 = ^Tr (d^Ud^) + ^Tr [ X U ] + U X ] 



where 



X 




u 



exp \i- 



3 

i=i 



Ti = 




V2 
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In the SU(2)l x SU(2)^ sector it is common to express quantities in the chiral limit without index 0, 
e.g., F and B. By this one means the SU(2)l x SU(2)ij chiral limit, i.e., m u = m d = but m s at its 
physical value. In the SU(3)l x SU(3)r sector the quantities Fq and Bq denote the chiral limit for all 
three quarks: m u = m d = m s = 0. Using the substitution U <-> IP, we see that £2 contains even powers 
of only: 

c 2 = cf + cf + ---. 

Since C 2 does not produce a vertex with three Goldstone bosons, there are no s-, u-, and t-channel 
pole diagrams, i.e., at D = 2, txtx scattering is entirely generated by a four-Goldstone-boson-interaction 
term. Expanding 

+ l F 2F 2 6F 3 + 24F 4 + '"' 
the interaction term Ct^ is identified as 



48F 2 



Tr([0, d^][(f>, + 2BTr(M(P i ) 
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Pb, b Pd, d 

v / 

Pa, a p c ,c 

Figure 3: Lowest-order Feynman diagram for 7T7t scattering. The vertex is derived from £2, denoted by 
2 in the interaction blob. 



We note that substituting F — > Fq, B — > Bq and the relevant expressions for and the quark- mass 
matrix M. the corresponding formula for SU(3)l x SU(3)r looks identical. Inserting (p = <fiiTi and 
working out the traces yields 

1 M 2 

£2 = iT^i&d^&d^j^j - <f> i <f>id fi <i> j d li <l>j) + 7^7™0t0t0j0j, 



where M 2 = 2Bm. The Feynman rule derived from for Cartesian isospin indices a, b, c, and d reads 
(see Fig. [3D 



M 



s-M 2 t-M 2 u-M 2 

OabOcd 1" ac 0M — h O a dVbc ^ 

i 



y , 2 (tabbed + 5 ac 5 M + 5 ad 5 bc ) (A a + A b + A c + A d ) , (120) 
where = p\ — M 2 and s, t, and u are the usual Mandelstam variables, 

s = (p a + Pb) 2 , t=(p a - p c ) 2 , u = {p a - p d ) 2 . 

In general, the T-matrix element for the scattering process TT a (p a ) + Tf b (p b ) — > n c (p c ) + Tr d (p d ) can be 
parameterized as 

T ab -> cd ( Pa , Pb ; Pc , Pd ) = 6 ab 6 cd A(s, t, u) + 5 ac 5 bd A(t, s, u) + 6 ad 5 bc A(u, t, s), (121) 

where the function A satisfies A(s, t, u) = A(s, u, t) |Weinberg, 1966| . Since the last line of the Feynman 
rule of Eq. ( 11201 ) disappears, if the external lines satisfy on-mass-shell conditions, at 0(q 2 ) the prediction 
for the function A is given by 

s - M 2 

A(s,t,u) = S —-^. (122) 

In Eq. (11221) we substituted F n for F and M n for M, because the difference is of 0(q A ) in T. Equation 
(11221) illustrates an important general property of Goldstone-boson interactions. If we consider the 
(theoretical) limit M 2 ,s,t,u — > 0, the T matrix vanishes, T — > 0. In other words, the strength of 
Goldstone-boson interactions vanishes in the zero-energy and mass limit. 

Usually, 7T7T scattering is discussed in terms of its isospin decomposition. Since the pions form an 
isospin triplet, the two isovectors of both the initial and final states may be coupled to / = 0,1,2. 
For m u = m d = rh the strong interactions are invariant under isospin transformations, implying that 
scattering matrix elements can be decomposed as 

(rX\T\I,h)=T J 5 u ,8 hI ,. (123) 
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For the case of tttt scattering the three isospin amplitudes are given in terms of the invariant amplitude 
A of Eq. fimi) by [Gasser and Leutwyler, 1984] 



T 
T 
T 



1=0 



1=1 



1=2 



3A(s, t, u) + A(t, u, s) + A(u, s, t), 
A(t, u, s) — A(u, s, t), 
A(t,u,s) + A(u,s,t). 



(124) 



For example, the physical tt + tt + scattering process is described by T I=2 . Other physical processes are 
obtained using the appropriate Clebsch-Gordan coefficients. 

Evaluating the T matrices at threshold, one obtains the s-wave 7T7r-scattering lengths 



->J=0| 



|thr 



32tto°, 



T 



I=2\ 



|thr 



327raQ, 



(125) 



where the subscript refers to s wave and the superscript to the isospin. (T 1 1 | thr vanishes because of 
Bose symmetry). The convention in ChPT differs by a factor (— M^) from the usual definition of a scat- 



tering length in the effective range expansion. The current-algebra prediction of Ref. Weinberg, 1966 
is identical with the lowest-order result obtained from Eq. (I122p . 



7M, 2 
32ttF2 



0.159, 



Ml 



16ttF 2 



-0.0454, 



(126) 



= 139.57 MeV. In order to 
Mthr = 0. Equations (11261) 



where we made use of the numerical values = 92 A MeV and M v = 
obtain the results of Eq. (11261) . use has been made of s t h r = 4M 2 and t t hr 
represent an absolute prediction of chiral symmetry. Once F n is known (from pion decay), the scat- 
tering lengths are predicted. The s-wave 7T7r-scattering lengths have been calculated at next-to-leading 
(NL) order Gasser and Leutwyler, 1984] and at next-to- next-to-leading order |Bijnens et al., 1 996 1 



Bij nens et al., 1997 . By matching the chiral representation of the scattering amplitude with a dis- 



persive representation |Roy 19 71] , |Ananthanarayan et al., 2001 , the predictions for the s-wave tttt 



scattering lengths are |Colangelo et al., 2000] , |Colangelo et al., 2001b] 



0.220 ±0.005, 



-0.0444 ±0.0010. 



(127) 



The empirical results for the s-wave 7T7r-scattering lengths have been obtained from various sources. In 



the K ei decay K 



it ' it e ' v e 



the connection with low-energy tttt scattering stems from a partial- 



wave analysis of the form factors relevant for the K ei decay in terms of tttt angular momentum eigen- 
states. In the low-energy regime the phases of these form factors are related by (a generalization of) 
Watson's theorem Watson, 1954 to the corresponding phases of I = s-wave and 1 = 1 p-wave elas- 
tic scattering |Colangelo et al., 2001a| . Using effective field theory techniques, isospin-breaking effects 
generated by real and virtual photons, and by the mass difference of the up and down quarks were 
discussed in Ref. |Colangelo et al., 2009| . 



Performing a combined analysis of the Geneva-Saclay data 



|Rosselet et al., 1977| , the BNL-E865 data |Pislak et al., 2001| , |Pislak et al., 2003| , and the NA48/2 



data |Batley et al., 2008 results in Colangelo et al., 2 009 



a u = 0.217 ± 0.008 cxp ± 0.006 th 



:i28) 



which is in excellent agreement with the prediction of Eq. (11271) . The — > TT ± TT + n reactions require an 
extrapolation to the pion pole to extract the tttt amplitude and are thus regarded to contain more model 
dependence, a° = 0.204 ± 0.014 (stat) ± 0.008 (syst) |Kermani et al., 1998| . The DIRAC Collaboration 
|Adeva et al., 2005 makes use of a lifetime measurement of pionium to extract |a§ — clq\ = 0.264^o;q2o- 



Finally, in the K — > Tr*^ ^ decay, isospin-symmetry breaking leads to a cusp structure ~ a® — a-i in 
the 7r°7r° invariant mass distribution near s^o^o rs 4M^+ Cabibbo, 2004 , Cabibbo and Isidori, 2 005 1. 
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Figure 4: Generic one-loop diagram. The black box denotes some unspecified vertex structure which is 
irrelevant for the discussion. 

Based on the model of |Cabibbo and Isidori, 2005] , the NA48/2 Collaboration extract a[J — a 2 , = 0.268 ± 
0.010 (stat) ±0.004 (syst) ±0.013 (ext). A more sophisticated analysis of the cusps in K — > 3n within an 
effective field theory framework can be found in Refs. |Colangelo et al., 2006| , |Bissegger et al., 2008 
and |Bissegger et al., 2009| . 



In particular, when analyzing the data of Ref. Pislak et al., 200T] in combination with the Roy equa- 



tions, an upper limit < 16 was obtained in Ref. Colangelo et al., 2001a| for the scale- independent 
low-energy coupling constant which is related to I3 of the SU(2)i x SU(2)# Lagrangian of Gasser and 
Leutwyler |Gasser and Leutwyler, 1984] . The great interest generated by this result is to be understood 
in the context of the pion mass at 0(q A ) 

Ml = M 2 — ^r^M 4 + 0(M 6 ), (129) 

where M 2 = 2mB. Recall that the constant B is related to the scalar quark condensate in the chiral 
limit and that a non-vanishing quark condensate is a sufficient criterion for spontaneous chiral sym- 
metry breakdown in QCD. If the expansion of M 2 in powers of the quark masses is dominated by 
the linear term in Eq. (11291) . the result is often referred to as the Gell-Mann-Oakes-Renner relation 
|Gell-Mann et al.~ 968|. If the terms of order m 2 were comparable or even larger than the linear terms, 
a different power counting or bookkeeping in ChPT would be required |Knecht et al., 1995] . The es- 
timate I/3 1 < 16 implies that the Gell-Mann-Oakes-Renner relation is indeed a decent starting point, 
because the contribution of the second term of Eq. (11291) to the pion mass is approximately given by 

T M 2 

4 = -0.054M^ for l 3 = 16, 



64tt 2 F 2 " 



i.e., more than 94 % of the pion mass must stem from the quark condensate Colangelo et al., 2001a . 



3.2.5 Primer to dimensional regularization 

If we want to use the Lagrangian of Eq. (11141) beyond the tree level, we will encounter ultraviolet diver- 
gences from loop integrals. For the regularization of the loop diagrams we will make use of dimensional 
regularization |'t Hooft and Veltman, 1972] , Leibbrandt, 1975] , |'t Hooft and Veltman, 1979 , because 



it preserves algebraic relations between the Green functions (Ward identities). We will illustrate the 
method by considering the following simple example, 

d 4 k i 

(130) 



(2vr) 4 A; 2 -M 2 + iO+ : 



which shows up in the generic diagram of Fig. HI Introducing a = yk 2 + M 2 > 0, we define f(ko) — 
{[ko + (a — *0 + )] [A; — (a — i0 + )]} _1 . In order to determine /f^ dk f(ko) as part of the calculation of /, 
we consider / in the complex ko plane and make use of Cauchy's theorem § c dzf(z) = for functions 
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rm(k )Aj 



-(a-iO+) 



Yl 



Y 3 t 



Y 4 V, 



a-iO+ 



Re(k ) 



Figure 5: Path of integration in the complex k plane. 

which are differentiable in every point inside the closed contour C. Choosing the path as shown in 
Fig. [5] and taking account of the fact that the quarter circles at infinity do not contribute, we obtain 
the so-called Wick rotation 



dtf(t)=-i dtf(it)=i dtf{it) 



(131) 



As an intermediate result, the integral of Eq. ( 1 1 3 f) reads 

d 4 / 1 



(2tt) 4 I 2 + M 2 ' 



where I 2 = ^1+^2 + ^3 + ^4 denotes a Euclidian scalar product in four dimensions. Performing the 
angular integration in four dimensions and introducing a cutoff A for the radial integration, the integral 
/ diverges quadratically for large values of / (ultraviolet divergence): 



1(A) 



dl 



A 2 M 2 , 



M 2 



8ir 2 Jo l 2 + M 2 (An) 2 (An) 2 \A 2 + M 2 / 

M 2 rl 



(132) 



where x = M 2 /A 2 — > as A — > 00. In dimensional regularization, we generalize the integral from 4 to 
n dimensions and introduce polar coordinates 

h = lcos(9i), I2 = I sin(6'i) cos^), h = I sin(^i) sin(6*2) cos(^3), 



l n -i = I sin(6 l i) sin( 



cos(6 l „_i), l n = Zsin(^i) sm(9 2 ) ■ ■ ■sin(9 n ^ 1 ) 



where < I, 9i G [0, vr] (i = 1, • • • , n — 2), and n _i G [0, 2n]. A general integral is then symbolically of 
the form 



d n l 



dll n ~ l / d9 n _ x / d0„_ 2 sin(0 n _ 2 ) 
Jo Jo 



d9i sin 



n~2i 



If the integrand does not depend on the angles, the angular integration can explicitly be carried out: 



7T 2 

d^ri 2 z r- 

r (l) 
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We define the integral for n dimensions (n integer) as 



I n (M 2 ^ 2 )=^- n J 



d n k i 
(2vr)" k 2 - M 2 + i0+ 



where the scale \x ('t Hooft parameter, renormalization scale) has been introduced so that the integral 
has the same dimension for arbitrary n. The integral formally reads 



In(M 2 ,/! 2 



4— nr 



7T 2 



dl- 



l 



n-l 



M 2 



4-n 



r (2) (2tt)» Jo l 2 + 

2^^- - — (M ^~ — — — i - 

r(f) (2tt)« 2 K 



r(i) 



4-n 



(4tt 



2\S-li 



r l- 



(133) 



Since T(z) is an analytic function in the complex plane except for poles of first order in 0, —1, —2, • • •, 
and a z = exp[ln(a)z], a G R + is an analytic function in <C, the right-hand side of Eq. (11331) can be 
thought of as a function of a complex variable n which is analytic in C except for poles of first order for 
n = 2, 4, 6, • • •. The analytic continuation for complex n reads 



I(M 2 ,fi 2 ,n) 



M 2 



'4:7111 



2\ 2 -2 



r i 



M 2 
16t? 



R + ln 



' M 2 ^ 
J*. 



where 



R 



n- 4 



[ln(4vr) + r(l) + l]. 



+ 0(n-4), (134) 



(135) 



The comparison between Eqs. f!134j) and (11321) illustrates the following general observations: in dimen- 
sional regularization power-law divergences are analytically continued to zero and logarithmic ultraviolet 
divergences of one-loop integrals show up as single poles in e = 4 — n. 



3.2.6 Power-counting scheme 

The Lagrangian £ c g of mesonic chiral perturbation theory is organized as a string of terms with an 
increasing number of derivatives and quark-mass terms, 

C eS = C 2 + U + A + • " • , (136) 

where the subscripts refer to the order in the momentum and quark-mass expansion. The index 2, for 
example, denotes either two derivatives or one quark-mass term. In terms of Feynman rules, derivatives 
generate four-momenta. A quark-mass term counts as two derivatives because of Eqs. (11051) - (11071) 
[M 2 ~ m q ) in combination with the on-shell condition p 2 = M 2 . We will generically count a small 
four-momentum — or the corresponding derivative — and a Goldstone-boson mass as of 0(q). The chiral 
orders in Eq. (11361) are all even [0(q 2k ), k > 1], because Lorentz indices of derivatives always have to 
be contracted and quark-mass terms count as 0(q 2 ). 

Besides the knowledge of the most general Lagrangian, we need a method which allows one to assess 
the importance of different renormalized diagrams contributing to a given process. For that purpose 
we analyze a given diagram under a simultaneous re-scaling of all external momenta, pi \— > tpi, and 
the light-quark masses, m q \— > t 2 m q (corresponds to M 2 i— > t 2 M 2 ). The chiral dimension D of a given 
diagram is defined as 

M(t Pi ,t 2 m g ,tfx) = t D M(pi,m g ,fx) = 0{q D ). (137) 
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D = 4-2-2-3 + 2-2 = 6. D = 4-2-2-3+2-1+4-1 = 8. D = 4-4-2-5 + 2-2 = 10. 

Figure 6: Application of the power-counting formula of Eq. (11381) in n = 4 dimensions. 

For small enough momenta (and masses) contributions with increasing D become less important. The 
chiral dimension is given by 

oo 

D = nN L - 2JV> + 2kN ^ (138) 



fc=i 

oo 



= 2 + (n-2)N L + Y.2(k-l)N 2k (139) 

k=l 

> 2 in 4 dimensions, 

where n is the number of space-time dimensions, Nl the number of independent loops, Ni the number 
of internal Goldstone boson lines, and N 2k the number of vertices from L 2k . Equation (11391) establishes 
a relation between the momentum and loop expansion, because at each chiral order, the maximum 
number of loops is bounded from above. In other words, we have a perturbative scheme in terms of 
external momenta and masses which are small compared to some scale [here: A x = AttFq — O (1 GeV)]. 
Examples of the application of the power-counting formula are shown in Fig. [61 

In order to prove the power-counting formula we start from the Feynman rules for evaluating the 
S-matrix element and investigate the behavior of the individual building blocks. Internal lines are 
described by a propagator in n dimensions which under re-scaling behaves as 

d n k i f d n k i k = tk> tn - 2 f d n k' 



(2n) n k 2 — M 2 + i0 + J (2ir) n t 2 (k 2 /t 2 — M 2 + i0 + ) J (2n) n k' 2 — M 2 + i0 + 

Vertices with 2k derivatives or k quark- mass terms re-scale as 

5 n (q)q 2k -> t 2k - n 5 n (q)q 2k , 

since p — > tp if q is an external momentum, and k = tk' if q is an internal momentum (see above). These 
are the rules to calculate S ~ 5 n (P)Ai. We need to add n to compensate for the overall momentum- 
conserving delta function. Applying these rules, the scaling behavior of the contribution to M. of a 
given diagram reads 

oo 

D = n + (n - 2)Ni + £ N 2k (2k - n). 

k=l 

The relation between the number of independent loops, the number of internal lines, and the total 
number of vertices N v = Y^kLi ^2k is given by Nl = Nj — (N v — 1). The product of N v momentum- 
conserving 8 functions contains overall momentum conservation. Therefore, one has Ny — 1 rather than 
Ny restrictions on the internal momenta. Applying 

oo 

-n £ N 2k = -nN v = n(N L — A/ — 1) 

k=l 
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| 2 for fc = 1. v ; | 4 for k = 1. 

Figure 7: The loop diagram is only suppressed if k m i n > 0. 

results in Eq. fl 13Sj) : 

oo 

D = nN L - 2Nj + J2 2kN 2k . 

k=i 

On the other hand, applying 

oo oo 

-n N 2k = -2 E N 2k + {n-2) (N L -Nj-1), 

k=l k=l 

results in Eq. ffT31?j) : 

oo 

5 = 2 + ^2(fc- l)N 2k + (n - 2)iV L . 
fe=i 

In particular, diagrams containing loops are suppressed due to the term 2Nl in four dimensions. 

Note that a minimal > is important. Otherwise, an infinite number of diagrams containing 
vertices from Cq would have to be summed (see Fig. [7]). This is for example the case when dealing with 
the nucleon-nucleon interaction. 



3.3 Next-to-leading order 



Already in 1967 it was shown by Weinberg Weinberg, 1967] that an effective Lagrangian is a con- 



venient tool for reproducing the results of current algebra in terms of tree-level calculations. In the 
purely mesonic sector, £ 2 of Eq. (I114p represents the corresponding Lagrangian. It was noted by Li 



and Pagels Li and Pagels, 1971| that a perturbation theory around a symmetry which is realized in 



the Nambu-Goldstone mode, in general, leads to observables which are non-analytic functions of the 
symmetry-breaking parameters, here the quark masses. In 1979 Weinberg initiated the application of 
an effective-field-theory program beyond the tree level allowing for a systematic calculation of correc- 
tions to the chiral limit [Weinberg, 1979] . When calculating one-loop graphs, using vertices from £ 2 , 
one generates ultraviolet divergences which in the framework of dimensional regularization appear as 
poles at space-time dimension n = 4. The loop diagrams are renormalized by absorbing the infinite 
parts into the redefinition of the fields and the parameters of the most general Lagrangian. Since C 2 
is not renormalizable in the traditional sense, the infinities cannot be absorbed by a renormalization 
of the coefficients F and B . However, to quote from Ref. [Weinberg, 1995] : "... the cancelation of 
ultraviolet divergences does not really depend on renormalizability; as long as we include every one of 
the infinite number of interactions allowed by symmetries, the so-called non-renormalizable theories are 
actually just as renormalizable as renormalizable theories." According to Weinberg's power counting 
of Eq. (I139p . one-loop graphs with vertices from £ 2 are of 0(q 4 ). The conclusion is that one needs to 
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construct the most general Lagrangian £4 and to adjust (renormalize) its parameters to cancel one- loop 
infinities. 

Beyond the quantum corrections to processes already described by £2, at next-to-leading order we 
will encounter another important feature, namely, the effective Wess-Zumino-Witten (WZW) action. 
The WZW action provides an effective description of the constraints due to the anomalous Ward 
identities. In general, anomalies arise if the symmetries of the Lagrangian at the classical level are not 
supported by the quantized theory after renormalization. 



3.3.1 The 0(q i ) Lagrangian of Gasser and Leutwyler 

The most general SU(3) L xSU(3) R -invariant Lagrangian at 0(q 4 ) is given by |Gasser and Leutwyler, 1985 



A = Li {Tr^C/fTW)*]} 2 + L 2 Tr [d^U(D u U) ] ] Tr [d»U(D v U) 



+L 3 Tr 
+L 5 Tr 
+L 7 



D^U(D»U) ] D V U(D V U) ] ] + L 4 Tr \D „U [D^] Tr ( X U ] + U X 



D,U(D»Uy( X U j + UxK 



Tr ( x rt + U X ] 



Tr ( X tft - Urf)\ +L 8 Tr(Ux j U X j + xU j X U j 



(140) 



The numerical values of the low-energy coupling constants Lj are not determined by chiral symmetry. 
In analogy to F and B of £2 they are parameters containing information on the underlying dynamics 
and should, in principle, be calculable in terms of the (remaining) parameters of QCD, namely, the 
heavy-quark masses and the QCD scale Aqcd- in practice, they parameterize our inability to solve 
the dynamics of QCD in the non-perturbative regime. So far they have either been fixed using empiri- 



cal input or theoretically using QCD-inspired models, meson-resonance saturation |Ecker et al., 1989a 
Pich, 2008| , and lattice QCD (see Ref. [Necco, (2009)] for a recent overview . 



By construction Eq. (11401) represents the most general Lagrangian at C(g 4 ), and it is thus possible 
to absorb the one-loop divergences by an appropriate renormalization of the coefficients Lj and Hf. 



r, 



32tt 2 



tR 



10), Hi 



1,2) 



(141) 



where R has already been defined in Eq. H135[) : 



n 



[ln(4vr) + r(l) + l], 



with n denoting the number of space-time dimensions and 7^ = — T'(l) being Euler's constant. The 
constants Tj and A$ are given in Table [TJ Except for L 3 and L7, the low-energy coupling constants Lj 
and the "contact terms" — i.e., pure external field terms — Hi and H 2 are required in the renormalization 
of the one-loop graphs. Since Hi and H 2 contain only external fields, they are of no physical relevance. 
The idea of renormalization consists of adjusting the parameters of the counter terms of the most 
general effective Lagrangian so that they cancel the divergences of (multi-) loop diagrams. In doing so, 
one still has the freedom of choosing a suitable renormalization condition. For example, in the minimal 
subtraction scheme (MS) one would fix the parameters of the counter term Lagrangian such that they 
would precisely absorb the contributions proportional to 2/(n — 4) in R, while the modified minimal 
subtraction scheme of ChPT (MS) would, in addition, cancel the term in the square brackets. 



38 



Coefficient 


Empirical Value 


r* 


L\ 


0.4 ±0.3 


3 

32 


L r 2 


1.35 ±0.3 


3 
16 


rr 
L 3 


-3.5 ± 1.1 





rr 
^4 


-0.3 ±0.5 


1 


rr 
^5 


1.4 ±0.5 


8 


rr 
-^6 


-0.2 ±0.3 


if 
144 


rr 
Li 7 


-0.4 ±0.2 





rr 
^8 


0.9 ±0.3 


5 

48 


rr 
L 9 


6.9 ±0.7 


1 
4 


rr 
^10 


-5.5 ±0.7 


1 
4 



Table 7: Renormalized low-energy coupling constants L\ in units of 10 3 at the scale \i 
|Bijnens et al., 1995| . Ai = -1/8, A 2 = 5/24. 



M p , see 



The renormalized coefficients L\ depend on the scale fi introduced by dimensional regularization 
[see Eq. (11341) ] and their values at two different scales Hi and /x 2 are related by 



L^2) = m^) + -^M-\- (142) 




We will see that the scale dependence of the coefficients and the finite part of the loop-diagrams 
compensate each other in such a way that physical observables are scale independent. 

A discussion of the two-flavor Lagrangian at 0(q A ) [Gasser and Leutwyler, 1984] can be found in 
Appendix D of Ref. Scherer, 2003| . For the construction of the 0(q 6 ) Lagrangian of even intrinsic 



parity, see Refs. Scherer and Fearing, 1995] , |Fearing and Scherer, 19 96] and |Bijnens et al., 1999 . For 



a status report on mesonic chiral perturbation theory beyond the one-loop level, we refer the reader to 
Ref. [Bijnens, 2007] . 



3.3.2 The effective Wess-Zumino-Witten action 

The Lagrangians discussed so far have a larger symmetry than QCD [Witten, 198 3]. For example, if we 
consider the case of "pure" QCD, i.e., no external fields except for the quark-mass term x = 2BqJ\4, £ 2 
and £4 contain interaction terms with an even number of Goldstone bosons only (even intrinsic parity). 
In other words, they cannot describe, e.g, K + K~ — > 7r + 7r~7r°. Analogously, £ 2 and £4 including a 
coupling to electromagnetic fields cannot describe 7r° — > 77. 

In order to overcome this shortcoming, Witten suggested to add the simplest term possible which 
breaks the symmetry of having only an even number of Goldstone bosons at the Lagrangian level. For 
the case of massless Goldstone bosons without any external fields the modified equation of motion reads 

9fl {-f Ud ^ U ^ + ^ upa Ud^U ] Ud v U ] Ud p U ] Ud a U ] = 0, (143) 

where A is a (purely imaginary) constant. For the purpose of writing down an action corresponding to 
Eq. (11431) . we extend the range of definition of the fields to a hypothetical fifth dimension, 

U(y) = exp ( ia^- ) , y* = (x M , a), i = 0, ■ • • , 4, < a < 1, (144) 
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where Minkowski space is defined as the surface of the five-dimensional space for a — 1. The action in 



in the absence of external fields (denoted by a superscript 0) is given by |Witten, 1 983 



= nS° wzw , S° wzw = -^- 2 f Q da j d 4 *e^ m Tr (lif • • <) , (145) 



where 



dll , in 



601234 = -6° 1234 = 1, Ut = U^, X 



dy i ' 48tt 2 ' 

A rather unusual and surprising feature of Eq. (I145P is that the action functional corresponding to the 
new term cannot be written as the four-dimensional integral of a Lagrangian expressed in terms of U 
and its derivatives. Wess and Zumino derived consistency or integrability relations which are satisfied by 
the anomalous Ward identities and then explicitly constructed a functional involving the pseudoscalar 
octet which satisfies the anomalous Ward identities |Wess and Zumino, 1971 . In particular, Wess and 



Zumino emphasized that their interaction Lagrangians cannot be obtained as part of a chiral invariant 
Lagrangian. Using topological arguments Witten showed that the constant n appearing in Eq. (11451) 
must be an integer. However, it was pointed out in Ref. |Bar and Wiese, 2001 that the traditional 



argument relating n with the number of colors N c is incomplete. Before discussing this argument, let 
us investigate the consequences of S^zw 

Expanding the SU(3) matrix U(y) in terms of the Goldstone boson fields, U(y) = 1 + ia<f)(x)/F + 
0(<j) 2 ), one obtains an infinite series of terms, each involving an odd number of Goldstone bosons, i.e., 
the WZW action 5' wzw is of odd intrinsic parity. For each individual term the a integration can be 
performed explicitly resulting in an ordinary action in terms of a four-dimensional integral of a local 
Lagrangian. For example, the term with the smallest number of Goldstone bosons reads 

^zw = 240^5 / d 4 xe^ CT Tr(0^0^09 p 09 ff 0), (146) 

which describes, e.g., K + K~ — > 7r + 7r~7r°. In particular, the WZW action without external fields involves 
at least five Goldstone bosons Wess and Zumino, 1971| . 



The connection to the number of colors N c is established by introducing a coupling to electromag- 



netism Wess and Zumino, 1971 , Witten, 1983 . In the presence of external fields there will be an 



additional term in the anomalous action, 

•Sano = S® no + S^ Q = n(Swzw + Swzw)' (147) 
given by (see, e.g., Ref. |Bijnens,T9 93|) 



with 



Swzw = J d 4 x e^Tr( Z^) (148) 



+iUd^lJ p U j r a - ilPd^TpUla + id p r u Ul p U ] r a - id p l v U ] r p Ul a 
-iUftvtfrpUl,, + iU^r v Ul P U ] r a - iU^U p l a + iU*r u r p r a 

+^A%rfd v T p Ul a - hj*Ud u l p U ] r a + ^U^rMdpL - ^U*Ul u U ] d p r a 
-Uffltfr p Ul a +UKU*Ul p tfr a + \ U^l v U L p l a - \ U*r v Ufr a 



+U^l u d p l a - U*r v d p r a + Utd v l P la - Uad„r p r a 

tfu^h + iU*U*U l p 



iU L M L v U% + ilAWUfr^ (149) 
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where we defined the abbreviations U p = Wd^U and = Ud^W. 

As a special case, let us consider a coupling to external electromagnetic fields by inserting 

where Q is the quark-charge matrix. The terms involving three and four electromagnetic four-potentials 
vanish upon contraction with the totally antisymmetric tensor & vpo ', because their contributions to Zp VpcT 
are symmetric in at least two indices, and we obtain 

2 

^wzw = -erU M J" + i^e^d v A p A a Tv[2Q 2 {Ud^ - U ] d p U) - QU^Qd^U + QUQd p U\ (150) 
We note that the current 

,fj,vpa 

JM = 4g~~2 ^(QdvUlfldpUlfldaUU* + QU ] d u UU ] d p UU ] d„U), e i23 = 1, (151) 

by itself is not gauge invariant and the additional terms of Eq. (j!50p are required to obtain a gauge- 
invariant action. The standard procedure of determining n is to investigate the interaction Lagrangian 
which is relevant to the decay 7r° — > 77 by expanding U = 1 + zdiag(7r°, — 7r°, 0)/F o + ■ ■ • . However, 
as pointed out by Bar and Wiese, when considering the electromagnetic interaction for an arbitrary 
number of colors one should replace the ordinary quark charge matrix by 

(I o\ f^r c + \ 
Q= o-l - jfc-i 

v° -u V 5^.-1 

The corresponding effective Lagrangian for 7r° — > 77 decay, 

n e 2 „„„„ „ „ 7T° 



results in the decay rate 



r o = f^ pa T t 

77 77 AL327T 2 ^ pu 



^ a 2 M\ n 2 „ , r f n 

r -''-=64^ivr 7 - 6eVx Gv; 



in good agreement with the experimental value (7.7 ± 0.6) eV for n = N c . However, the result is no 



indication for N c = 3 Bar and Wiese, 2001 . The conclusion from their analysis is that one should 
rather consider three-flavor processes such as 77 — > 7r + 7r~7 or — > i^7r to test the expected N c 
dependence in a low-energy reaction. For example, the Lagrangian relevant to the decay rj — > 7r + 7r~7 is 
given by 

16-Tl 

£rpr+7r- 7 = — r- „A Qu - Qd)^ vp0 r A p d u r]dp7i + d a 7v' , 

where the quark-charge difference Q u — Qd = 1 is independent of N c . However, by investigating the 
corresponding rj and rj' decays up to next-to-leading order in the framework of the combined 1/N C and 
chiral expansions, Borasoy and Lipartia have concluded that the number of colors cannot be determined 
from these decays due to the importance of sub-leading terms which are needed to account for the 
experimental decay widths and photon spectra |Borasoy and Lipartia, 2 005 . 



For a discussion of the 0(q ) Lagrangian of odd intrinsic parity see Refs. |Ebertshauser et al., 2002 
and |Bijnens et al., 2002| . 



41 



Figure 8: Self-energy diagrams at 0(q 4 ). Vertices derived from C 2n are denoted by 2n in the interaction 
blobs. 



= + ®— - + + 

Figure 9: Unrenormalized propagator as the sum of irreducible self-energy diagrams. Hatched and cross- 
hatched "vertices" denote one-particle-reducible and one-particle-irreducible contributions, respectively. 



3.3.3 Masses of the Goldstone bosons 

A discussion of the masses at 0(q A ) is one of the simplest applications of chiral perturbation theory 
beyond the tree level. For that purpose let us consider £ 2 + £4 for QCD with finite quark masses but in 
the absence of external fields. We restrict ourselves to the limit of isospin symmetry, i.e., m u = = m. 
In order to determine the masses we calculate the self energies £Qo 2 ) of the Goldstone bosons. 
Let 

A ^M = p ,-M| 2 + . i0+ - * = < 152 > 
denote the Feynman propagator containing the lowest-order masses, 

M 2 2 = 2B m, M| i2 = B (m + m s ), M 2 2 = ^B (m + 2m s ) . 

(The subscript 2 refers to chiral order 2.) The proper self-energy insertions, — i'E^p 2 ), consist of one- 
particle- irreducible diagrams only, i.e., diagrams which do not fall apart into two separate pieces when 
cutting an arbitrary internal line. At chiral order D = 4, the contributions to — zE^ 4 (p 2 ) are those 
shown in Fig. [HJ In general, the full (unrenormalized) propagator may be summed using a geometric 
series (see Fig. |H]): 

i 

p 2 -M2 2 -£ (p 2 )+*O+' 
The physical mass, including the interaction, is defined as the pole of Eq. f)153j) . 

MJ - M 2 2 - S (M 2 ) ± 0, (154) 

where the precision of the determination of M| depends on the precision of the calculation of S^,. 

For our particular application with exactly two external meson lines, the relevant interaction La- 
grangians can be written as 

Ant = C-i + 1 (155) 

where 

^ = ^F2{Tr([0,^0]0^0) + J B o Tr(^0 4 )}, (156) 



— Mi 



i0- 



+ 



(153) 
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Figure 10: Contribution of the pion loops to the 7T° self energy. 



-a K K + K~ - b K d,K + d^K~ - a K K°K° - b K d„K°d»K c 



-- [a v r] + b^rid^rij . 



The constants au and bs are given by 



a w = 


QAB 2 r 


b n = 


16£? 
F 2 


a K = 


3252 r 

F 2 i 


b K = 


Fi 


a v = 


UBl r 

3F 2 L 


b v = 


16£? 



(2m + m s )mL G + m 2 L 5 
[(2m + m s )L4 + mLs] 



(2m + m s )(rh + m s )L 6 + -(m + m s ) 2 L s 



(2m + m s )L 4 + -(m + m s )L 5 



(2m + m s )(m + 2m s )L 6 + 2(m — m s ) 2 L 7 + (m 2 + 2m 2 )L i 



(2m + m s )L 4 + ~(m + 2m s )L 5 



At C?(g 4 ) the self energies are of the form 



£<m(p 2 



A 6 + 



(157) 



(158) 



(159) 



where the constants A§ and receive a tree-level contribution from £ 4 and a one-loop contribution 
with a vertex from £ 2 (see Fig. [8]). For the tree-level contribution of £ 4 this is easily seen, because the 
Lagrangians of Eq. (11571) contain either exactly two derivatives of the fields or no derivatives at all. For 
example, the contact contribution for the 77 reads 



-i^:(p 2 ) = -i(a v + b vP 2 ). 

For the one-loop contribution the argument is as follows. The Lagrangian contains either two 
derivatives or no derivatives at all which, symbolically, can be written as 00<90<90 and 4 , respectively. 
The first term results in M 2 or p 2 , depending on whether the <fi or the dcj) are contracted with the 
external fields. The "mixed" situation vanishes upon integration. The second term, 4 , does not 
generate a momentum dependence. 

As a specific example, we evaluate the pion-loop contribution to the tc° self energy (see Fig. [10]) by 
applying the Feynman rule of Eq. f|120p for a = c = 3, p a = p c = p, b = d = j, and Pb = Pd = k: 
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d 4 k 



(2vr) 4 3F- 



■[-Ap 2 -4k 2 + 5il£ 



7T,2J 



k 2 -M 2 



7T,2 



«0 H 



where 1/2 is a symmetry factor. Since the integral diverges, we consider its extension to n dimensions. 
In addition to the loop-integral of Eq. (11341) , we need 

d n k k 2 f d n k k 2 -M 2 + M 2 



(2ir) n k 2 - M 2 + i0+ 



(2tt)™ k 2 - M 2 + iO- 



where we have added 



-M 2 + M 2 in the numerator. We make use of 



d n k 
(2vr)" 







in dimensional regularization which is "shown" as follows. Consider the (more general) integral 

d n k(k 2 ) p , (160) 



substitute k = Xk' (A > 0), and relabel k' — k 

= \ n+2p I d n k{k 2 ^ p 



Since A > is arbitrary and, for fixed p, the result is to hold for arbitrary n, Eq. ( 11601) is set to zero 
in dimensional regularization. We emphasize that the vanishing of Eq. ( 1 1 6 f) has the character of a 
prescription. The integral does not depend on any scale and its analytic continuation is ill defined in 
the sense that there is no dimension n where it is meaningful. It is ultraviolet divergent for n + 2p > 
and infrared divergent for n + 2p < 0. 
We then obtain 

d n k k 2 

M 2 I(M 2 ^ 2 ,n), 



with I(M 2 ,fi 2 ,n) of Eq. 



(2n) n k 2 -M 2 + i0+ 
The pion-loop contribution to the 7r° self energy is thus 

' -Ap 2 + M 2 )I(M 2 fx 2 ,n), 



6F 2 



which is indeed of the type discussed in Eq. ( 11591) and diverges as n — > 4. 

After analyzing all loop contributions and combining them with the tree-level contributions of Eqs. 
(11581) . the constants A$ and of Eq. (I159P are given by 



Btt 

A k 
B k 

An 



K2 

Fo 2 



6 



i« 2 



1 I(M 2 t2 ) - \l{Mlt) +32[(2m + m s )B L 6 + mB L 8 ] 



6 



2I(M 2 



1 I(Mk,2 

F 2 + 3 F? 



loop contribution 

165(5 r/n - 

[{2m + m s )L 4 + mL 5 



tree-level contribution 











F 2 



F 2 ' 
1I(M- 
I' 

F 2 

Fo 2 



^« 2 ) - \l(M 2 t2 ) - \l{M 2 Kt2 ) + 32 



(2m + m s )B L 6 + -(m + m s )B L 8 



17,2) 



+ 



1 1(M 2 . 



2) , ll(Mh) 



F 2 



+ 



+ 



K2 

F n 2 



F 2 



16- 



B 
'f 2 



{2m + m s )L i + -{m + m s )L 5 



2 ) 

r/,2) 



MM 2 ) + -I{M 2 



[16M 2 2 L 8 + 32(2m + m s )B L 6 ] + 



128 B 2 (m 



m. 



F 2 



-(3L 7 + L 8 



B r , 



I{M\ 



F 2 



16 , , M. 

—2 (2m + m s )B L 4 - 8- 



n,2 



F 2 U > 



(161) 
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where, for simplicity, we have suppressed the dependence on the scale \i and the number of dimensions 
n in the integrals I(M 2 ,fi 2 ,n) [see Eq. (I134p ]. Both the integrals / and the bare coefficients Li (with 
the exception of L 7 ) have l/(n — 4) poles and finite pieces. In particular, the coefficients and 
are not finite as n — > 4 showing that they do not correspond to observables. 

The masses at 0(q 4 ) are determined by solving Eq. (1154ft with the predictions of Eq. (11591) for the 
self energies, 

Ml 



M^ + A^ + B^Ml 



from which we obtain 



Ml 



Mj 2 + A, 
1-B* 



MUl + B^+A^ + Oiq 6 



because A$ = 0(q 4 ) and {B^, Af| 2 } = 0(q 2 ). Expressing the bare coefficients L, in Eq. (11611) in terms 
of the renormalized coefficients by using Eq. (1141 j) . the results for the masses of the Goldstone bosons 
at C(g 4 ) read Gasser and Leutwyler, 1985] 



ML 

7T,4 



M 2 KA 



M, 



q A 



<2 1 + 



16 



Ml 2 (M^ 



32tt 2 F 2 



96tt 2 F 2 



+ -=2 [(2m + m a )B (2L r 6 - L\) + mB (2L r 8 - L r 5 ) 



(162) 



M 2 \l + 



M l ln (Ki 



48vr 2 F 2 



16 
+ F? 







(2m + m s )B (2L r 6 - L\) + -(m + m s )B (2L r 8 - L 



(163) 
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16tt 2 F 2 
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2An 2 F 2 
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+M\ 



2 (2m + m a )S (2L5 - L4) + 8-^(2L^ - L5) 
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96vr 2 F 2 
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32tt 2 F 2 
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'tt,2 



M|, 2 
48tt 2 F 2 



In 



' M 



K,2 



128 5 2 (m 



9 



-(3L^ + Lg). 



(164) 



First of all, we note that the expressions for the masses are finite. The infinite parts of the coefficients 
Li of the Lagrangian of Gasser and Leutwyler exactly cancel the divergent terms resulting from the 
integrals. This is the reason why the bare coefficients Li must be infinite. Furthermore, at 0(q i ) the 
masses of the Goldstone bosons vanish, if the quark masses are sent to zero. This is, of course, what we 
had expected from QCD in the chiral limit but it is comforting to see that the self interaction in £. 2 ( m 
the absence of quark masses) does not generate Goldstone boson masses at higher order. At 0(q 4 ), the 
squared Goldstone boson masses contain terms which are analytic in the quark masses, namely, of the 
form m 2 multiplied by the renormalized low-energy coupling constants L[. However, there are also non- 



analytic terms of the type m ln(m g ) — so-called chiral logarithms- 



Such a behavior is an illustration of the mechanism found by Li and Pagels Li and Pagels, 1971 



which do not involve new parameters. 

who 



noticed that a perturbation theory around a symmetry which is realized in the Nambu- Goldstone mode 
results in both analytic as well as non-analytic expressions in the perturbation. Finally, the scale 
dependence of the renormalized coefficients L\ of Eq. (11411) is by construction such that it cancels the 
scale dependence of the chiral logarithms. Thus, physical observables do not depend on the scale fi. It 
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is straightforward to verify this statement by differentiating Eqs. (11621) - (j!64p with respect to fi and 
by making use of 

dLm i\_ 

d/i 167r 2 /x' 

where the Tj are given in Table [7J 

3.3.4 Electromagnetic polarizabilities of the pion 

Another strong constraint provided by chiral symmetry is the connection between the electromagnetic 
polarizabilities of the charged pion and the radiative pion beta decay In the framework of classical 
electrodynamics, the electric and magnetic polarizabilities a and j3 describe the response of a system to 
a static, uniform, external electric and magnetic field in terms of induced electric and magnetic dipole 
moments. In principle, empirical information on the pion polarizabilities can be obtained from the 
differential cross section of low-energy Compton scattering on a charged pion, 



da fu'V e 2 \ e 2 1 + z 



uuj r 



(a + [3) n+ (l + z) 1 + (a - (3) w +(l - z)' 



+ 



where z = q ■ q' and uj'/uj = [1 + u(l — z)/M n }. The forward and backward differential cross sections 
are sensitive to (a + (3) n + and (a — respectively. 

Within the framework of the partially conserved axial- vector (PCAC) hypothesis and current algebra 
the electromagnetic polarizabilities of the charged pion are related to the radiative charged-pion beta 
decay 7r + — > e + z/ e 7 |Terent'ev, 1973| . The result obtained using ChPT at leading non-trivial order 



(0(q 4 )) |Bijnens and Cornet, 1 988 is equivalent to the original PCAC result, 

e 2 1 I 



'Aix (4vrF 7r ) 2 M^ 6 



where /a = (^6 — ^5) is a linear combination of scale-independent parameters of the two-flavor 0(q A ) La- 
grangian |Gasser and Leutwyler, 1 984|. At 0(q A ) this difference is related to the ratio 7 = Fa/ ' Fy 



of the pion axial-vector form factor Fa and the vector form factor Fy of radiative pion beta de- 
cay |Gasser and Leutwyler, 1984| , 7 = /a/6. Once this ratio is known, chiral symmetry makes an abso- 
lute prediction for the polarizabilities. This situation is similar to the s-wave 7T7r-scattering lengths of 
Eq. (11261) which are predicted once F n is known. Using the most recent determination 7 = 0.443 ±0.015 
by the PIBETA Collaboration |Frlez et al., "20041 ( assuming Fy = 0.0259 obtained from the conserved 
vector current hypothesis) results in the 0(q A ) prediction a n + = (2.64 ± 0.09) x 10 _4 fm 3 , where the 
estimate of the error is only the one due to the error of 7 and does not include effects from higher orders 
in the quark-mass expansion. 

Corrections to the leading-order PCAC result have been calculated at 0(q 6 ) and turn out to be 
rather small |Biirgi, 1996| , Gasser et al., 2006] . Using updated values for the LECs, the predictions of 



|Gasser et al., 2006 ] are 

(a + p)^ = 0.16 x 10~ 4 fm 3 , (165) 
(£*-/?)*+ = (5.7 ± 1.0) x KT 4 fm 3 . (166) 

The corresponding corrections to the 0(q 4 ) result indicate a similar rate of convergence as for the tttt- 
scattering lengths |Gasser and Leutwyler, 1984] , |Bijnens et al., 1996 . The error for (a + [3) n + is of the 



order 0.1 x 10" 4 fm , mostly from the dependence on the scale at which the 0(q 6 ) low-energy coupling 
constants are estimated by resonance saturation. 

As there is no stable pion target, empirical information about the pion polarizabilities is not easy to 
obtain. For that purpose, one has to consider reactions which contain the Compton scattering amplitude 
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Figure 11: The reaction 7p — > "fK + n contains Compton scattering on a pion as a sub diagram in the t 
channel, where t = (jp n — p p ) 2 - 




Figure 12: Differential cross section averaged over 537 MeV < E y < 817 MeV and 1.5 M% < s\ < 5M^. 
Solid line: model 1; dashed line: model 2; dotted line: fit to experimental data. 

as a building block, such as, e.g., the Primakoff effect in high-energy pion-nucleus bremsstrahlung, 
n~Z — > 7r~Zj, radiative pion photoproduction on the nucleon, jp — > •yTt+n, and pion pair production 
in e + e~ scattering, e + e~ — ► e + e~7r + 7r~. Unfortunately, at present, the experimental situation looks 
rather contradictory (see Refs. Ahrens et al., 2005] , |Gasser et al., 2006] for recent reviews of the data 



and further references to the experiments). 

The potential of studying the influence of the pion polar izabilities on radiative pion photoproduction 
from the proton was extensively studied in |Drechsel and Fil'kov, 1994 . In terms of Feynman diagrams, 



the reaction yp — > 77r + n contains real Compton scattering on a charged pion as a pion pole diagram (see 
Fig. ITT]) ■ In the recent experiment on 7p — > / ~f7T + n at the Mainz Microtron MAMI Ahrens et al., 20 05 1, 



the cross section was obtained in the kinematic region 537 MeV < < 817 MeV, 140° < 023 < 180°. 
Figure [12] shows the experimental data, averaged over the full photon beam energy interval and over 
the squared pion-photon center-of-mass energy s\ from 1.5 Ml to 5 Ml ClS db function of the squared 
pion momentum transfer t in units of M%. For such small values of si, the differential cross section is 
expected to be insensitive to the pion polarizabilities. Also shown are two model calculations: model 1 
(solid curve) is a simple Born approximation using the pseudoscalar pion-nucleon interaction including 
the anomalous magnetic moments of the nucleon; model 2 (dashed curve) consists of pole terms without 
the anomalous magnetic moments but including contributions from the resonances A(1232), P n (1440), 
-D 13 (1520) and S f n(1535). The dotted curve is a fit to the experimental data. 

The kinematic region where the polarizability contribution is biggest is given by 5M^ < s\ < 15M^ 
and — YIMl < t < — 2M%. Figure [TBI shows the cross section as a function of the beam energy integrated 
over s\ and t in this second region. The dashed and solid lines (dashed-dotted and dotted lines) refer 
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Figure 13: The cross section of the process *yp —> ^7r + n integrated over s\ and t in the region where 
the contribution of the pion polarizability is biggest and the difference between the predictions of the 
theoretical models under consideration does not exceed 3 %. The dashed and dashed-dotted lines are 
predictions of model 1 and the solid and dotted lines of model 2 for (a — (3)^+ = and (a — f3) n + = 
14 x 10~ 4 fm 3 , respectively. 



to models 1 and 2, respectively, each with (a — (3)^+ = ((a — f3) n + = 14 x 10~ 4 fm 3 ). By comparing 
the experimental data of the 12 points with the predictions of the models, the corresponding values of 
(a — (3) n + for each data point have been determined in combination with the corresponding statistical 
and systematic errors. The result extracted from the combined analysis of the 12 data points reads 
Ahrens et al., 2005| 



(a - (3) w+ = (11.6 ± 1.5 stat ± 3.0 syst ± 0.5 mod ) x 10" 4 fm 3 (167) 

and has to be compared with the ChPT result of (5.7 ± 1.0) x 10 _4 fm 3 , which deviates by 2 standard 
deviations from the experimental result. 

Clearly, the model-dependent input to the result of Eq. fl 1 6 7|) deserves further study. In partic- 
ular, the model error was estimated by comparing the analysis with two specific models. In Ref. 
Kao et al., 2007| radiative pion photoproduction was studied in the framework of heavy-baryon chi- 



ral perturbation theory at the one-loop level. Unfortunately, the kinematical conditions of the MAMI 
experiment were not explicitly considered. It was argued that the extraction of pion polarizabilities 
is, in principle, possible and that the main uncertainty in the extraction arises from the effect of two 
structures of the 0(q 3 ) Lagrangian. 

The Primakoff method was used at Serpukhov with the result Antipov et al., 1983| 



(a - (3) n+ = (13.6 ± 2.8 stat ± 2.4 syst ) x 10~ 4 fm 3 , (168) 

in agreement with the value from MAMI. Recently, also the COMPASS Collaboration at CERN has 
investigated this reaction, and the data analysis is underway |Guskov, 2008] . Unfortunately, the third 
method based on the reactions e + e~ — > 77 — > 7r + 7r~, has led to even more contradictory results (see 
Ref. |Gasser et al., 2006| ). 

Also on the theoretical side there has been a long-standing problem. The application of disper- 
sion sum rules as performed in FiFkov and Kashevarow, 1999| , |Fil'kov and Kashevarow, 2006| yields 



a — P) n + = (13.0^9) x 10~ 4 fm which provides an even more pronounced discrepancy with the pre- 
dictions of chiral perturbation theory than the MAMI result |Gasser et al., 20 06 1. These dispersion 
relations are based on specific forms for the absorptive part of the Compton amplitudes. In Ref. 
|Pasquini et al., 2008 , the analytic properties of these forms have been examined and the strong en- 



hancement of intermediate-meson contributions was shown to be connected with spurious singularities. 
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It was shown that the results of dispersion theory and effective field theory are not in conflict, once the 
basic requirements of dispersion relations are taken into account. 

Clearly, the electromagnetic polarizabilities of the charged pion remain one of the challenging topics 
of hadronic physics in the low-energy domain. Chiral symmetry provides a strong constraint in terms 
of radiative pion beta decay and mesonic chiral perturbation theory makes a firm prediction beyond 
the current algebra result at the two-loop level. Both the experimental determination as well as the 
theoretical extraction from experiment require further efforts. 



4 Baryonic chiral perturbation theory 
4.1 Lagrangian 

So far we have considered the purely mesonic sector involving the interaction of Goldstone bosons with 
each other and with the external fields. Now we want to describe matrix elements with a single baryon 
in the initial and final states. 



4.1.1 Transformation properties of the fields 



Our aim is the most general description of the interaction of baryons with Goldstone bosons and external 
fields at low energies. For that purpose we not only need to specify the transformation behavior of the 
Goldstone bosons and external fields but also of the remaining dynamical fields entering the Lagrangian. 
Our discussion follows Refs. |Georgi, 1984| , Gasser et al., 1988] . Consider the nucleon doublet and the 
octet of | + baryons (see Fig. [2]), 
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n 



h V2 



V 



Lyo i J_ a 

^0 



p 

n 



(169) 
(170) 



Each entry of \P and B, respectively, is a complex, four-component Dirac field. In contradistinction 
to the case of the Goldstone-boson matrix of Eq. (1551) . we have B ^ B< . The representation of the 
isospin group SU(2)y and the flavor group SU(3)y on {9} and {-£>}, respectively, is given by 



^ i-> V^f, V e SU(2) V , 
B i-> VBV\ V e SU(3) y , 



(171) 
(172) 



i.e., \& transforms under the fundamental representation of SU(2) and B transforms under the adjoint 
representation of SU(3). Starting from Eqs. (1171 j) and (j!72j) we will discuss realizations of SU(2) L x 
SU(2) fl and SU(3) L x SU(3) fl on } and {5}, respectively. 

Let us begin with G = SU(2) L x SU(2) fl . Recall that the transformation of Eq. (j55|) . 

U i — > RUL\ 

defines a nonlinear realization of G on {U}. Introducing u 2 = U, we define the SU(2)-valued function 
K{L,R,U) by 



u 



u 



\fmW = RuK-\L,R,U), i.e., K(L, R, U) = u'^Ru = VRUtf l RVU. (173) 
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The transformation 

defines an operation of the group G on the set {(U,^f)}. This is true, because (a) the identity of G 
leaves any pair (U, \&) invariant and (b) the transformation satisfies the homomorphism property 

/ s / J u \ , J r 2 ul\ \ ( r x r 2 ul\l\ \ 

where we made use of 

K{L±, R x , R 2 UL\)K(L 2 , R 2 , U) = K{{L X L 2 ), (R 1 R 2 ), U). 

Note that for a general group element g = (L, R) G G the transformation behavior of ^ depends on 
U. The exception to this rule is the case of an isospin transformation R = L = V , where, because of 
U' = u' 2 = VuVWuV^ = Vu 2 V^ = VUV\ one has v! = VWL Comparing with Eq. ffT73l . we obtain 
K~ l (V, V, U) = or K(V,V,U) = V. This is consistent with our starting point that ^ transforms 
linearly as an isospin doublet under the isospin subgroup H = SU(2) y of G = SU(2) L x 811(2)^. Recall 
that the symmetry of the vacuum determines the multiplet structure of the spectrum |Coleman, 1 966 . 
For G = SU(3) L x SU{3) R one uses 

^ : (b)^( U b>) = ( k(l, r, u)bk1(l, r,u))> (175) 

where K is defined completely analogously to Eq. (1 1 73 j) after inserting the corresponding SU(3) matrices. 

The generalization to other multiplets is straightforward. One first specifies the transformation 
behavior under the subgroup H in terms of V and V^. In order to find the transformation behavior 
under G one simply replaces V — ► K and — > . 



4.1.2 Baryonic effective Lagrangian at lowest order 

Given the dynamical fields of Eqs. f)174p and (j!75j) and their transformation properties, we will now 
discuss the most general effective baryonic Lagrangian at lowest order. We will start with the effective 
ttN Lagrangian C^jj which we demand to have a local SU(2)i x SU(2) R x U(l) v symmetry. The 
transformation behavior of the external fields is given in Eq. (f49j) . whereas U and the nucleon doublet 
transform as 

( U(x) \ ( V R (x)U(x)Vl(x) \ 

1 J V ex P [-tQ(x)}K[V L (x),V R (x),U(x)}<f(x) I ' 



(176) 



The local character of the transformation implies that we need to introduce a covariant derivative 
with the usual property that it transforms in the same way as 

U M * = (^ + r^ (177) 

where the so-called connection is given by 



r M - 2 



u^dfj, - ir^u + u{d^ - il^)u ] . (178) 
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Since K not only depends on Vl and Vr but also on U, the covariant derivative contains besides the 
external fields also u and v) and their derivatives. At 0(q) there exists another Hermitian building 
block, the so-called vielbein, 



u^d^-ir^u-u^-il^ , (179) 



which under parity transforms as an axial vector, \— > — w M , and under SU(2) L x SU(2) fl x U(l)y, 
transforms as u M i— > Ku^K^ . 

The structure of the most general effective 7riV Lagrangian describing processes with a single nucleon 
in the initial and final states is of the type tyO^, where O is an operator acting in Dirac and flavor space, 
transforming under SU(2) L x SU(2) R x U(l)y as KOK^ . The Lagrangian must be a Hermitian Lorentz 
scalar which is even under the discrete symmetries C, P, and T. The most general such Lagrangian 
with the smallest number of derivatives is given by |Gasser et al., 1 988 1 

dH = ^(iip-m + |V75?v) *■ (180) 

It contains two parameters not determined by chiral symmetry: the chiral limit m of the nucleon 
mass m^v and the chiral limit g A of the axial- vector coupling constant g&. The physical value of is 
determined from neutron beta decay and is given by gA = 1-2695 ± 0.0029. The overall normalization 
of the Lagrangian is chosen such that in the case of no external fields and no pion fields it reduces to 
that of a free nucleon of mass m. 

Similarly as in the mesonic case, the Lagrangian C^ N has predictive power once the two pa- 
rameters have been identified. For example, a tree-level calculation of pion-nucleon scattering pro- 



duces the famous Weinberg- Tomozawa relation for the s-wave 7riV-scattering lengths Weinberg, 1966 



Tomozawa, 1966 



a 1 



87T(l+/i)F2 
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where 1 = 1/2 or 7 = 3/2 refers to the total isospin of the ttN system. As in ttit scattering, the s-wave 
7riV-scattering lengths vanish in the chiral limit, i.e., Goldstone bosons interact "weakly" with other 
hadrons in the zero-energy and mass limit. 

Since the nucleon mass does not vanish in the chiral limit, the zeroth component d° of the 
partial derivative acting on the nucleon field does not produce a "small" quantity. This results in 
new features of the chiral power counting in the baryonic sector. The counting of the external fields 
as well as of covariant derivatives acting on the mesonic fields remains the same as in mesonic chiral 
perturbation theory. On the other hand, the counting of bilinears tyTty is probably easiest understood 
by investigating the matrix elements of positive-energy plane-wave solutions to the free Dirac equation 
in the Dirac representation: 

^ + \t,x)=exp{-ip-x)jE + m N ( g| V (182) 



where \ denotes a two-component Pauli spinor and p M = (E,p) with E = yp 2 + m%. In the low- 
energy limit, i.e., for non-relativistic kinematics, the lower (small) component is suppressed as \p |/mjv 
in comparison with the upper (large) component. For the analysis of the bi-linears it is convenient to 
divide the 16 Dirac matrices into even and odd ones, £ = {l,7o,757tj ^ij} and O = {75, 7570, 7i, cno} 
|Foldy and Wouthuysen, 1950| , respectively, where odd matrices couple large and small components but 
not large with large, whereas even matrices do the opposite. Finally, id^ acting on the nucleon solution 
produces p M which we write symbolically as p M = (mjv,0) + (E — mjsr,p), where we count the second 
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term as 0(q), i.e., as a small quantity. We are now in the position to summarize the chiral counting 



scheme for the (new) elements of baryon chiral perturbation theory |Krause, 1990 

* = O(q ), = O(q ), (lip - m)tt = 0(q), 



1,Tm:T57 m .^ = 0{q% 75 = 0(q), (183) 

where the order given is the minimal one. For example, 7^ has both an O(q ) piece, 70, as well as an 0(q) 
piece, 7j. Note that because of the additional spin degree of freedom the baryonic effective Lagrangian 
contains both odd and even chiral orders. A rigorous non-relativistic reduction may be achieved in 
the framework of the Foldy-Wouthuysen method |Foldy and Wouthuysen, 1950| , Fearing et al., 1994 



or the heavy-baryon approach | Jenkins and Manohar, 1991| , |Bernard et al., 19 92a 



The construction of the SU(3) L x 811(3)^ Lagrangian proceeds similarly except for the fact that the 
baryon fields are contained in the 3x3 matrix of Eq. ( 11701) transforming as KBK^ . As in the mesonic 
sector, the building blocks are written as products transforming as K ■ ■ ■ with a trace taken at the 
end. The lowest-order Lagrangian reads |Georgi, 1984 , Krause, 1990 



'-MB — 11 



B [ip - M ) B\ - -Tr (f^K, B}) - -Tr {B^j 5 [u^, B}) , (184) 



where Mo denotes the mass of the baryon octet in the chiral limit. The covariant derivative of B is 
defined as 

D fM B = d fM B + [T^,B], (185) 

with of Eq. (j!78p [for SU(3) L x SU(3) fl ]. The constants D and F may be determined by fitting the 
semi-leptonic decays B — > B' + e~ + v e at tree level [Borasoy, 1999| : 

D = 0.80, F = 0.50. (186) 
4.2 Renormalization and power counting 

In the following discussion we will restrict ourselves to the two-flavor case. The effective Lagrangian 
relevant to the one-nucleon sector consists of the sum of the purely mesonic and ttN Lagrangians, 
respectively, 

£ eff = C n + C nN = £f + 4 4 ) + • • • + 4i + 4 2 i + • • • , (187) 

which are organized in a derivative and quark-mass expansion. Tree-level calculations involving the sum 
4 2 ^ + 4iV reproduce the current algebra results. The higher-order Lagrangians of the ttN sector can 
be found in |Gasser et al., 1988] , Ecker and Mojzis, 1996 , Fettes et al., 2000 . When studying higher 



orders in perturbation theory in terms of loop corrections one encounters ultraviolet divergences. As a 
preliminary step, the loop integrals are regularized, typically by means of dimensional regularization. In 
the process of renormalization the counter terms are adjusted such that they absorb all the ultraviolet 
divergences occurring in the calculation of loop diagrams. This will be possible, because we include in 



the Lagrangian all of the infinite number of interactions allowed by symmetries Weinberg, 1995 . At 
the end the regularization is removed by taking the limit n — > 4. Moreover, when renormalizing, we 
still have the freedom of choosing a renormalization condition. As we will see, the power counting is 
intimately connected with choosing a suitable renormalization condition. 

4.2.1 The generation of counter terms 

Before discussing the power-counting problem and its solution, let us briefly recall the principles of 
the renormalization procedure which will then allow us to set up a consistent power counting. At 
the beginning, the Lagrangian is written down in terms of bare, i.e., unrenormalized parameters and 
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fields. In order to illustrate the procedure let us discuss C^ N of Eq. (11801) and consider the free part in 
combination with the ttN interaction term with the smallest number of pion fields, 



(i) 



where the subscript B denotes bare quantities. The renormalization is performed by expressing all the 
bare parameters and bare fields of the effective Lagrangian in terms of renormalized quantities (see, 
e.g., Refs. Collins, 1984| , Weinberg, 1995 for details). Introducing the renormalized fields through 



B 



J %B 



189) 



we express the field redefinition constants and \JZ$ and the bare quantities in terms of renormalized 
parameters: 



Zq, = l + 6Z v (m,g A ,gi,v), 

Z^ = 1 + SZj (m, g A , g h v) , 

m B = m(v) +5m(m,g A ,gi,is) , 

Sab = Sa(") + Sg A (m,g A ,gi,u) 



(190) 



where gi, i = 1, • • • , oo, collectively denote all the renormalized parameters which correspond to bare 
parameters gi B of the full effective Lagrangian of Eq. (1187p . The parameter v indicates the dependence 
on the choice of the renormalization condition. We emphasize that the usual choice m[y) = m, where 
m is the nucleon pole mass in the chiral limit, is only one among an infinite number of possibilities. 
Substituting Eqs. fl]ggj) and fTM?]) into Eq. (jlggj) . we obtain 



(i) 

■kN 



C 



basic 



+ ^ct + 



(191) 



with the so-called basic and counter-term Lagrangians, respectively, 



(192) 
(193) 



where we introduced the abbreviations 



5{m} 



bZ^ra + ZySm, 

sz^^K + zJ^ 

r \ r b 



F 



% + 



Sa, 



In Eq. (11921) . m, g A , and F denote the chiral limit of the physical nucleon mass, the axial- vector 
coupling constant, and the pion-decay constant, respectively. Expanding the counter-term Lagrangian 
of Eq. (11931) in powers of the renormalized coupling constants generates an infinite series. By adjusting 
the expansion coefficients suitably, the individual terms are responsible for the subtractions of loop 
diagrams. For example, the divergences occurring in dimensionally regularized one-loop calculations 
involving vertices of £2 are absorbed in the renormalization of the bare coefficients Lj [see Eq. (I14ip . 
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Figure 14: One-loop contributions to the nucleon self energy. The number 1 in the interaction blobs 
refers to 4jv- 



4.2.2 Power counting for renormalized diagrams 

In the following, whenever we speak of renormalized diagrams, we refer to diagrams which have been 
calculated with a basic Lagrangian and to which the contribution of the counter-term Lagrangian has 
been added. Counter-term contributions are typically denoted by a cross. One also says that the 
diagram has been subtracted, i.e., the unwanted contribution has been removed with the understanding 
that this can be achieved by a suitable choice for the coefficient of the counter-term Lagrangian. In this 
context we will adjust the finite pieces of the renormalized couplings such that renormalized diagrams 
satisfy the following power counting: a loop integration in n dimensions counts as q n , pion and fermion 
propagators count as q~ 2 and respectively, vertices derived from and 4*zv count as q 2k and 
q h , respectively. Here, q collectively stands for a small quantity such as the pion mass, small external 
four-momenta of the pion, and small external three-momenta of the nucleon. The power counting does 
not uniquely fix the renormalization scheme, i.e., there are different renormalization schemes leading to 
the above specified power counting. 



4.2.3 The power-counting problem 

In the mesonic sector, the combination of dimensional regularization and the modified minimal sub- 
traction scheme MS [see Eq. (11411) ] leads to a straightforward correspondence between the chiral and 
loop expansions. By studying the one-loop contributions of Fig. [TJ] to the nucleon self energy, we will 
see that this correspondence, at first sight, seems to be lost in the baryonic sector. 

In the following we will calculate the mass of the nucleon up to and including 0(q 3 ). As in the 
case of the Goldstone bosons, the physical mass is defined through the pole of the full propagator, but 
here at p = m^. In terms of the nucleon self energy S(|>) we will solve the equation 

m N — m — S(m7v) = 0, (194) 

where m denotes the nucleon mass in the chiral limit. 

According to the power counting specified above, we need to calculate the two types of one-loop 
contributions shown in Fig. [TH together with the corresponding counter-term contribution and a tree- 
level contribution. After renormalization, we would like to have the orders D = n-1 — 2-1 — 1-1 + 2-1 = 
n — 1 for the first loop diagram and n-1 — 2-l + l- l= n — 1 for the second loop diagram. 

The basic interaction Lagrangian obtained from expanding C^ N up to and including two pion fields 
reads 

42 = —^-fiM^ - i*/?xai. f*. 



The corresponding Feynman rules are given by 
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AF 2 



+ f J^abc^c 



(196) 



(2) 

From the next-to- leading-order ttN Lagrangian C n jy we only need one term, namely, 



resulting in the constant tree-level contribution 



En ree = -4 Cl M 2 . 



(197) 



:i98) 



(3) 

Moreover, there are no tree-level contributions from the Lagrangian C^. The second diagram of Fig 
is zero, because the contraction e aac = in the Feynman rule of Eq. (11961) vanishes. In dimensional 
regularization, the first diagram of Fig. [TH generates the contribution 



•3g^ .4-71 



'AF 2 



d n k 



— m 



(27r) n [(p - k) 2 - m 2 + i0+}(k 2 — M 2 + i0+) ' 



(199) 



Using {7^,7i/} = 2p M1/ , the numerator of the integrand is written as 



+ m)k 2 + (p 2 — m 2 



yielding the intermediate result 

3gi 



£ loop (f>) 



+ m)^~ n i 



AF 2 { 

-{f + m)M 2 ^~ n i J 

+ (p 2 - m 2 )^- n i 
d n k 



d n k 



(p — k) 2 — 



1 



(27t)" (p - A;) 2 - m 2 + i0+ 
d n k 1 



(2?r) n [(p - A;) 2 - m 2 + iO+](A; 2 — M 2 + t0 + ) 
d n k $ 
(2yr) n [(p - k) 2 -m 2 + i0+}(k 2 — M 2 + i0+) 

* 1 



(27r) n A; 2 — M 2 + i0 + 



(200) 



The last term in Eq. (I200p vanishes since the integrand is odd in k. We use the following convention 
for scalar loop integrals, 



(Pi 



?1, 



d n fc 



1 



(2tt)" [(fc + pi) 2 - m 2 + t0+] ■■■[(k + qi ) 2 — M 2 + z0+] 
The vector integral in the third line of Eq. fl200j) is determined using the ansatz 

d n A; k u 



(201) 



(2tt)« [(p - A;) 2 - m 2 + iO+](A; 2 — M 2 + i0+) 
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(202) 




Figure 15: Renormalized one- loop self-energy diagram. 



Multiplying Eq. H202[) by p^, one obtains for C, 



(j — 

2p 2 



In~I*+ (P 2 ~m 2 + M 2 )I NlT (-p, 0) 



(203) 



In terms of the above convention for the scalar loop integrals the loop contribution to the nucleon self 
energy reads 



£ loop (f0 



M 

4F2 



+ m)I N + {f + m)M 2 I Nw (-p, 0) 



-(p 2 — m 2 ) 



2p 2 



IN 



h + (P 2 



m 2 + M 2 )I Nn ( 



-p,0) 



(204) 



The explicit expressions for the integrals are given by 

M 2 



In*(p, 0) 



167T 2 
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167T 2 



R + \n 
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IN 



R + \n 



R + In — - 1 + 



167T 2 

p 2 -m 2 - M 2 

p2 



m 



, (M\ 2mM ' 
In — +— Ffl 

V m / 
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where is given in Eq. (11351) . Q is defined as 



n 



P 



m 



M 2 



2mM 



and 



V^ 2 - Hn (-ft 



- VH 2 ^ 

Vl — ft 2 arccos(— f2), 
V^ 2 - 1 In ffi + Vft 2 - 1 



?7T 



Vn 2 



i. 



ft < -i, 
-i < ft < i, 
i < ft. 



Because of the terms proportional to R, the result for the self energy contains divergences as n — > 4, so 
it has to be renormalized. The counter-term Lagrangian must produce structures which precisely cancel 
the divergences, because otherwise the result for the nucleon mass will not be finite. For convenience, 
we choose the renormalization parameter \i — m. 

In the modified minimal subtraction scheme MS all the contributions proportional to R are canceled 
by corresponding contributions generated by the counter-term Lagrangian of Eq. (11931) . but also by 
counter-term Lagrangians resulting from higher-order terms of Eq. (1187!) . Operationally this means 
that we simply drop all terms proportional to R and indicate the renormalized coupling constants 
by a subscript r. Again, this is possible, because we include in the Lagrangian all of the infinite 



number of interactions allowed by symmetries Weinberg, 1995 . The renormalized diagram is depicted 
in Fig. [151 where the cross generically denotes counter-term contributions. The MS-renormalized self- 
energy contribution then reads 



4F 2 



+ m )M 2 I Nn (-p,0)-(p 2 



m 



2p 2 



(p 2_ m 2 + M 2 )I r N ^_ P]0) _ K 



(206) 
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where the superscript r on the integrals means that the terms proportional to R have been dropped. 
Writing m — 2m + (i> — w) and comparing the first term of Eq. (12041) with Eq. (12061) . we note that 

among other terms, the MS renormalization involves (even in the chiral limit) an infinite renormalization 



yielding the relation between the bare and the renormalized mass |Gasser et al., 1988 

3g|r 
32vr 2 F 2 



m B = m+ -^%—m 3 R + 



Using 

j U-p,o) = - t ^ + ---, 

we see that the MS-renormalized self energy produces a contribution of 0(q 2 ) which is in conflict with 
the power counting assigned above. For a long time this was interpreted as the absence of a systematic 
power counting in the relativistic formulation of ChPT. 
We can now solve Eq. (I194j) for the nucleon mass, 



m N 



m + Y% r ee (m N ) + E^mjv) = m - 4c lr M 2 + E|° op (m w ). (207) 



We have for the difference — m = 0(q 2 ). Since our calculation is only valid up to 0(q 3 ), it is 
sufficient to determine S|° op (mAr) to that order. In fact, using arccos (—ft) = f + • • • , the expansion of 

1 / TTM \ 

= ^ + — + ••• , (208) 



167T 2 V m 



from which we obtain for the nucleon mass in the MS scheme Gasser et al., 1988 



m N = m- Ac lr M 2 + - ^lr^! (209) 

32vr 2 F 2 32vr 2 F 2 v ; 

The solution to the power-counting problem is the observation that the term violating the power count- 
ing, namely, the third on the right-hand side of Eq. (12091) . is analytic in the quark mass and can thus 
be absorbed in counter terms. In addition to the MS scheme we have to perform an additional finite 
renormalization. For that purpose we rewrite 

c lr = c 1 + 6c u 5 Cl = ^^r 2 +... (210) 



3> 



in Eq. (1209!) which then gives the final result for the nucleon mass at 0(q 

3p- 2 A/f 3 

mN = m -A Cl M 2 -^— 2 . (211) 

To summarize, we have shown that the validity of a power-counting scheme is intimately connected 
with a suitable renormalization condition. In the case of the nucleon mass, the MS scheme alone 
does not suffice to bring about a consistent power counting. We will shortly outline two methods, the 



infrared renormalization Becher and Leutwyler, 1999 and the extended on-mass-shell renormalization 



Fuchs et al., 2003a , which both produce a systematic power counting in a manifestly Lorentz-invariant 



framework. 

4.3 Solutions to the power-counting problem 
4.3.1 Heavy-baryon approach 

The first solution to the power-counting problem was provided by the heavy-baryon formulation of ChPT 



Jenkins and Manohar, 1991 , |Bernard et al., 1992a| . The basic idea consists in dividing an external 
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p p+q p' 

Figure 16: s-channel pole diagram of 7rN scattering. 



nucleoli four-momentum into a large piece close to on-shell kinematics and a soft residual contribution: 
p = mv + k p , v 2 = 1, v° > 1 [often v M = (1, 0, 0, 0)]. The relativistic nucleon field is expressed in terms 
of velocity-dependent fields, 

with 

M v = e +"-(l + H v = e +™-(l - i>)^i. 

Using the equation of motion for 7i v , one can eliminate 7i v and obtain a Lagrangian for M v which, to 



lowest order, reads Bernard et al., 1992a 



3? N =K{iv-D + g A S v ■ u)Af v + 0(l/m), S£ = \i^v v 



The result of the heavy-baryon reduction is a 1/m expansion of the Lagrangian similar to a Foldy- 
Wouthuysen expansion |Foldy and Wouthuysen, 1950 . In higher orders in the chiral expansion, the ex- 



pressions due to 1/m corrections of the Lagrangian become increasingly complicated |Fettes et al., 2 000 



Moreover — and what is more important — the approach sometimes generates problems regarding 
analyticity which can be illustrated by considering the example of pion-nucleon scattering |Becher, 2 002 



The invariant amplitudes describing the scattering amplitude develop poles for s = m 2 N and u = m 2 N . 
For example, the singularity due to the nucleon pole in the s channel (see Fig. [TBI) is understood in 
terms of the relativistic propagator 

! ' (212) 



(p + q) 2 - m 2 N 2p ■ q + M 2 ' 

which, of course, has a pole at 2p-q = —M 2 or, equivalently, s = m 2 N . Analogously, a second pole results 
from the u channel at u — m 2 N . Although both poles are not in the physical region of pion-nucleon 
scattering, analyticity of the invariant amplitudes requires these poles to be present in the amplitudes. 
Let us compare the situation with a heavy-baryon type of expansion, where, for simplicity, we choose 
as the four- velocity p^ = mjsrv 1 *, 

1 11 \ \ ( Ml \ , , 

m^ = ——A 1 - ^ : : +•••)■ (213) 



2p ■ q + 2m N v . n -\- Mjl. 2m N v ■ q \ 2m N v ■ q 

Z.TY1. tit * 



Clearly, to any finite order the heavy-baryon expansion produces poles at v ■ q = instead of a simple 
pole at v-q — —M 2 /(2m N ) and will thus not generate the (nucleon) pole structures of the invariant am- 
plitudes unless an infinite number of diagrams is summed. For a comprehensive overview of calculations 
performed in the heavy-baryon framework the reader is referred to Ref. |Bernard et al., 19 95 1 . 
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4.3.2 Master integral 

We have seen that the modified minimal subtraction scheme MS does not produce the desired power 
counting. We will discuss the power-counting problem in terms of the dimensionally regularized one-loop 
integral 

d n k 1 



H(p 2 ,m 2 ,M 2 ;n) = -i f 



(27r) n [(k - p) 2 - m 2 + i0+](k 2 — M 2 + i0+) 
. r d n k 1 

(2tt)" [k 2 - 2p ■ k + {p 2 - m 2 ) + i0+] {k 2 - M 2 + i0+) ' ^ ^ 

We are interested in nucleon four-momenta close to the mass-shell condition, p 2 m m 2 , counting p 2 — m 2 
as 0(q) and M 2 as 0(q 2 ). In order to conform with Ref. Becher and Leutwyler, 1999| , we have omitted 
the factor /i 4_n and have reversed the overall sign in comparison with our previous definition of I^ n . Let 
us turn to the discussion of H(p 2 , m 2 , M 2 ; n). To that end, we make use of the Feynman parametrization 

1 r 1 1 

dz l 7TT, \v> (215) 



ab Jo [az + b(l — z)} 2 

with a = (k — p) 2 — m 2 + i0 + and b = k 2 — M 2 + i0 + , interchange the order of integrations, and perform 
the shift k — > k + zp, to obtain 

d n k 1 



H(p 2 , m 2 , M 2 ; n) = -i j dz J 



where 

Making use of 



(2tt)™ [k 2 - A(z) + i0+} 2 ' 
A(z) = z 2 p 2 - z(p 2 -m 2 + M 2 ) + M 2 . 

d n k (k 2 y i(_)p-9r(p + f)r(g-p-f) 



(2nr(k 2 -A)« (4vr)f r (f ) T(q) 

with p = and q = 2, we find 

H(p 2 ,m 2 ,M 2 ;n) = -^rT (2-^-) [' dz[A{z) - z0+]§- 2 . (216) 
(Air) 2 \ 2/ Jo 

The relevant properties can nicely be displayed at the threshold p 2 hr = (m + M) 2 , where A(z) = 
[z(m + M) — M] 2 is particularly simple. The small imaginary part can be dropped in this case, because 
A(z) is never negative. Splitting the integration interval into [0,z ] and [z , 1] with z Q = M/(m + M), 
we have, for n > 3, 

dz[A(zp~ 2 = / dz[M - z(m + M)} n ~ A + / dz[z{m + M) — M] n ~ 4 

JO Jzq 

(n-3)(m + M) v y ' 
yielding, through analytic continuation, for arbitrary n 



H((m + M) 2 , m 2 , M 2 ; n) = 7 ^ 27 : + ^—-^ . (217) 

VV ^' ' ' ; (4vr)2(n-3) V^ + ^f m + MJ K J 

The first term, proportional to M n_3 , is defined as the so-called infrared singular part /. Since M — > 
implies p 2 hr — > m 2 this term is singular for n < 3. The second term, proportional to m n ~ 3 , is defined as 
the infrared regular part R. Note that for non-integer n the infrared singular part contains non-integer 
powers of M, while an expansion of the regular part always contains non-negative integer powers of M 
only. 
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4.3.3 Infrared regularization 

Let us now turn to a formal definition of the infrared singular and regular parts for arbitrary p 2 



Becher and Leutwyler, 1999] which makes use of the Feynman parametrization of Eq. (12161) . Introduc- 



ing the dimensionless variables 

M ~, N ^ p 2 -m 2 -M 2 in/ nN 

a = — = 0(q), n = P = O(q ), 218 

m 2mM 

we rewrite A(z) as 

A{z) = m 2 [z 2 - 2attz{l - z) + a 2 (l - z) 2 } = m 2 C(z), 
so that H is now given by 

H{p 2 , m 2 , M 2 - n) = n(m; n) dz[C(z) - i0 + ]%~ 2 , (219) 

Jo 

where 



r (2 — - 

K (m;n)= ) x „ 2y m n ~ 4 . (220) 

(47r) 2 

The infrared singularity originates from small values of z, where the function C(z) goes to zero as 
M — > 0. In order to isolate the divergent part one scales the integration variable z = ax so that the 
upper limit z = 1 in Eq. (12191) corresponds to x = 1/a oo as M -> 0. An integral / having the same 
infrared singularity as H is then defined which is identical to H except that the upper limit is replaced 
by oo: 

poo poo 

I = K(m; n) / dz[C{z) - i0 + ] ~ 2 = fc(m; n)a n ~ 3 / dx[D(x) - i0 + ] " 2 , (221) 
J Jo 

where 

D(x) = 1 - 2Qx + x 2 + 2ax(Qx - 1) + a 2 x 2 . 
(The pion mass M is not sent to zero.) Accordingly, the regular part of H is defined as 

poo 

R=-K{m-n) J dz[C{z) -i0 + }%- 2 , (222) 

so that 

H = I + R. (223) 

Let us verify that the definitions of Eqs. (12211) and (I222p indeed reproduce the behavior of Eq. (I217p . 
To that end we make use of Q t hr = 1, yielding 

ithr = K{m;n)a n ~ 3 J dx {[(1 + a)x - l] 2 - i0 + \ 2 , (224) 



which converges for n < 3. In order to continue the integral to n > 3, write Becher and Leutwyler, 1999 



and make use of a partial integration 

J dx {{{1 + a)x - l] 2 -z0+} 2 = 

~ (l + a)x-l fr/ _ , .„+-if-2l°° 1 f°° . fr „ , 1 §-2 



-- {[(1 + a)x - l] 2 - z0+} 2 "^jf d2; {[( 1 + «) a; - 1 ] 2 -^ + } 



_> 



(l + a)(n-4) 
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For n < 3, the first expression vanishes at the upper limit and, at the lower limit, yields l/[(l+a)(n— 4)]. 
Bringing the second expression to the left-hand side, we may then continue the integral analytically as 



I d, {[(! + <,), -!]'-*)*}■ = , n _ 3)(1 + g) . P25) 



so that we obtain for 7 t hr 



'thr 



K(m;n)a n - 3 - - = - — \, ' , 226 

V ' ; (n-3(l + a) 4vr 2 n-3 m + M' v ' 



which agrees with the infrared singular part I of Eq. (121 7p . 

The threshold value of the regular part of Eq. (12221) is obtained by analytic continuation from n < 3 
to n > 3: 



r(2-f) roc 

R thT = Tr^ir 1 / + M) — M n 

(47r) 2 Jl 



n— 4 



r( 2 -f) 



(4?r)3 (n-3)(m + M) 



r^3 r ( 2 ~f) m-3 

(47r)f (n- 3) m + M' 1 J 

which is indeed the regular part R of Eq. (12171) . 

What distinguishes I from R is that, for non-integer values of n, the chiral expansion of I gives 
rise to non-integer powers of 0(q), whereas the regular part R may be expanded in an ordinary Taylor 
series. For the threshold integral, this can nicely be seen by expanding 1^ and R t hr in the pion mass 
counting as 0(q). On the other hand, it is the regular part which does not satisfy the counting rules. 
The basic idea of the infrared renormalization consists of replacing the general integral H of Eq. (12161) 
by its infrared singular part J, defined in Eq. (I22ip . and dropping the regular part R, defined in Eq. 
In the low-energy region H and / have the same analytic properties whereas the contribution of 



.R, which is of the type of an infinite series in the momenta, can be included by adjusting the coefficients 
of the most general effective Lagrangian. This is the infrared renormalization condition. 

As discussed in detail in Ref. Becher and Leutwyler, 1999] , the method can be generalized to an 



arbitrary one-loop graph (see also Ref. Semke and Lutz, 2006| ). It is first argued that tensor integrals 



involving an expression of the type k^ 1 ■ ■ ■ k^" 2 in the numerator may always be reduced to scalar loop 
integrals of the form 

A n k 1 1 



-i 



(27r)™ a x • ■ -a m bi • • -b n 

where a« = (qi + k) 2 — M 2 + i0 + and &j = (pi — k) 2 — m 2 + i0 + are inverse meson and nucleon propagators, 
respectively. Here, the qi refer to four- momenta of 0(q) and the pi are four-momenta which are not 
far off the nucleon mass shell, i.e., p 2 = m 2 + 0(q). Using the Feynman parametrization, all pion 
propagators and all nucleon propagators are separately combined, and the result is written in such a 
way that it is obtained by applying (m — 1) and (n — 1) partial derivatives with respect to M 2 and m 2 , 
respectively, to a master formula. A simple illustration is given by 

V L . /",,: 1 



0102 Jo [axz + a 2 (l — z)} 2 DM 2 Jo a\z + a 2 (l — z) ' 

where = (qi + k) 2 - M 2 +i0 + . Of course, the expressions become more complicated for larger numbers 
of propagators. The relevant property of the above procedure is that the result of combining the meson 
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propagators is of the type If A with A = (k + q) 2 — M 2 + i0 + , where q is a linear combination of the m 
momenta g^, with an analogous expression 1/B for the nucleon propagators. Finally, the expression 



_.rd n k 1 
~J (2ir) n AB 

may then be treated in complete analogy to H of Eq. (I214p . i.e., the denominators are combined 
as in Eq. (I215p . and the infrared singular and regular pieces are identified by writing dz ■ ■ ■ = 



4.3.4 Extended on-mass-shell scheme 

In the following, we will concentrate on yet another solution which has been motivated in Ref. 



Gcg elia and Japaridze, 1999 and has been worked out in detail in Ref. |Fuchs et al., 2003a . The cen- 



tral idea of the extended on-mass-shell (EOMS) scheme consists of performing additional subtractions 
beyond the MS scheme such that renormalized diagrams satisfy the power counting. Terms violating 
the power counting are analytic in small quantities and can thus be absorbed in a renormalization of 
counter terms. 

In order to illustrate the approach, let us consider a simplified version of the integral H, namely its 
value in the chiral limit, 



H (p 2 , m 2 , 0; n) = —i J 



d n k 



1 



(2n) n [{k - p) 2 -m 2 + i0+](k 2 + i0+) 



where 



A 



2 2 

p — m 

m 2 



0{q) 



is a small quantity. Applying the power-counting rules of Sec. I4.2.2[ we want the renormalized integral 
to be of order D = n — 1 — 2 = n — 3. Introducing C(z, A) = z 2 — A z(l — z) — i0 + , we obtain 



H(p 2 ,m 2 ,0;n) 



dz [C(z,A)Y 



(228) 



where K,(m; n) is given in Eq. (12201) . For the purpose of evaluating the integral of Eq. f!228j) we write 
£ dz [C(z, A)]f- 2 = (-A)f- 2 jT 1 dz z%- 2 (l 



1 + A \ 2- 2 



A 



and apply Eqs. 15.3.1 and 15.3.4 of Ref. Abramowitz and Stegun, 1972] to obtain 



H(p 2 ,m 2 ,0;n) = K{m; w) ^V^- F ( 1, n nn .V 



2 ' 2 ' m 2 r 



(229) 



where F(a,b;c;z) is the hypergeometric function |Abramowitz and Stegun, 1972| . In order to discuss 
the power-counting properties of H (in the chiral limit) in terms of A, we make use of Eq. 15.3.6 of 
Ref. |Abramowitz and Stegun, 1972] to re-write Eq. (I229p as 



H(p 2 ,m 2 ,0;n) 



m 



n—A 



(4tt) 



r(2-§) / 

— ^ -P-F 1, 2 - -: 4 - n: - 



n 



A 



+(-A)- 3 r g - l) r(3 - n)F Q - 1, n - 2; n - 2; - a) 



(230) 
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Making use of 



, s oh a(a + l)b(b + l)z 2 

F(a, b;c;z) = l + -z+ 1 ; \ J - + 
c cc+1 2 



(231) 



for | z | < 1 and the fact that A counts as a small quantity of 0(q), we immediately see that the first 
term of Eq. ( 12301) contains a contribution which does not satisfy the above power counting, i.e., which 
is not proportional to 0(q) as n — > 4. Using the expansion of Eq. (12311) together with T(l +x) — xT(x) 



we obtain, as n 



H 



m 



n— 4 



(47T) 



r(2-f 



n 



P 



P 



P 



P 



l_i_ l n l_i_ + l_^_ k 1-1- 



77T 



(232) 



where • ■ ■ refers to terms which are at least of 0(q 3 ) or 0(n — 4). Note that we count a term of the 
type — Aln(— A) as 0(q). If we subtract 



m 



n— 4 



r(2-f 



from Eq. 



(47r) 2 n — 3 
we obtain as the renormalized integral 



(233) 



H R (p 2 ,m 2 ,0;n) 



m 



n— 4 



( 47r )n/2 




V 77T , 



+ ■ 



(234) 



The subtracted term of Eq. (I233j) is local in the external momentum p, i.e., it is a polynomial in p 2 and 
can thus be obtained by a finite number of counter terms in the most general effective Lagrangian. 
We have seen in Eq. (12301) that the one-loop integral is of the type 



#~F(n,A) + A n - 3 G(n,A), 

where F and G are hypergeometric functions and are analytic in A for any n. The observation central 
for the setting up of a systematic method is the fact that the part proportional to F can be obtained 
by first expanding the integrand in small quantities and then performing the integration for each term 
|Gegelia et al., 1994] (see Sec. 14.3.51 for an illustration of the general method). We now apply a conven- 
tional renormalization prescription which allows us to identify those terms which we subtract from a 
given integral without explicitly calculating the integral beforehand. In essence we work with a modified 
integrand which is obtained from the original integrand by subtracting a suitable number of counter 
terms. The meaning of suitable in the present context will be explained in a moment. To that end we 
consider the series 



E 

1=0 



(p 2 



m 



2\l 



/! 



1 







1 



2p 2P "dp fl ) [k 2 -2k-p+ {p 2 - m 2 ) + i0+]{k 2 + 10 + ) 



(k 2 

+ (p 2 — m 2 



2k ■ p + i0+)(k 2 + i0+) 
1 1 



1 



2m 2 (k 2 - 2k ■ p + i0+) 2 2m 2 (k 2 - 2k ■ p + i0+)(k 2 + i0+) 
1 



(k 2 -2k-p + i0 + ) 2 (k 2 + i0~ 



+ ■ ■ 



(235) 



where [• • ■] p 2 =m 2 means that we consider the coefficients of (p 2 — m 2 ) 1 only for four-momenta p^ 1 satisfying 
the on- mass-shell condition. Although the coefficients still depend on the direction of after integra- 
tion of this series with respect to the loop momentum k and evaluation of the resulting coefficients for 
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p 2 = m 2 , the integrated series is a function of p 2 only. In fact, as was shown in Ref. Gegelia et al., 1994) , 
the integrated series exactly reproduces the first term of Eq. (12301) . At this point we stress that 



—i 



d n k 



(2tt)" (k 2 - 2k ■ p + i0+)(k 2 + i0 H 



and 



r d n k 1 
J (2vr)« (k 2 - 2k ■ p + p 2 - m 2 + i0+)(k 2 + ifr 



are not the same for n < 3. 

The formal definition of the EOMS renormalization scheme is then as follows: we subtract from the 
integrand of H(p 2 ,m 2 ,0;n) those terms of the series of Eq. (I235j) which violate the power counting. 
These terms are always analytic in the small parameter and do not contain infrared singularities. In 
the above example we only need to subtract the first term. All the higher-order terms contain infrared 
singularities. For example, the last term of the second coefficient would generate a behavior k 3 /k 4 of 
the integrand for n = 4. The integral of the first term of Eq. (12351) is given by Eq. (I233j) . and we end 
up with Eq. (I234[) for the renormalized integral: 



Hr = H — H, 



subtr 



0(q r 



Since the subtraction point is p 2 



m 



the renormalization condition is denoted "extended on-mass- 



shell" (EOMS) scheme in analogy with the on-mass-shell renormalization scheme in renormalizable 
theories. In the general case including the pion mass, one would consider the series 



1 



(A; 2 - 2k ■ p + i0+) (k 2 + i0+) 



+M Z 



ik 2 - 2k ■ p + i0+) ik 2 + to+y 



+ (p 2 — m 2 ) 
+ ••• 



1 



2m 2 [k 2 -2k-p + i0- 



+ 



instead of Eq. ( 1235h . However, it would still be only the contribution resulting from the first term that 
were to be subtracted. 

Within the EOMS framework it is straightforward to obtain a consistent power counting in mani- 



festly Lorentz-invariant baryon ChPT including, e.g., vector mesons Fuchs et al., 2003b| or the A(1232) 
Hacker et al., 2005] as explicit degrees of freedom. Moreover, the infrared regularization 



resonance 



of Becher and Leutwyler can be reformulated in a form analogous to the EOMS renormalization 
scheme and can thus be applied straightforwardly to multi-loop diagrams with an arbitrary number of 
particles with arbitrary masses |Schindler et al., 2004a| (see also Refs. |Lehmann and Prezeau, 2 002 1, 
|Bruns and Meifiner, 2005] , Bruns and Meifiner, 2008] ). The application of both infrared and extended 
on-mass-shell renormalization schemes to multi-loop diagrams was explicitly demonstrated by means of 
a two-loop self-energy diagram Sc hindler et al., 2004b| . 



4.3.5 Dimensional counting analysis 



In this section we provide an illustration of the dimensional counting analysis |Gegelia et al., 1994 



in 



terms of a specific example. To that end let us consider the one- loop integral of Eq. (I214p . 



r d n k 1 1 

H(p 2 ,m 2 } M 2 ;n) = -i J j—y n - _ - - - ~ 2 _ — - ^ + k2 _ M2 + %Q+ - (236) 
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One would like to know how the integral behaves for small values of M and/or p 2 — m 2 as a function 
of n. If we consider, for fixed p 2 ^ m 2 , the limit M —> 0, the integral H can be represented as 

H(p 2 , m 2 , M 2 ; n) = J2 M^F^p 2 , m 2 , M 2 ; n), (237) 

i 

where the functions Fi are analytic in M 2 and are obtained as follows. First, one re- writes the integration 
variable as k = M ai k, where a.i is an arbitrary non-negative real number. Next, one isolates the overall 
factor of so that the remaining integrand can be expanded in positive powers of M 2 and interchanges 
the integration and summation. The resulting series represents the expansion of Fi(p 2 ,m 2 ,M 2 ;n) 
in powers of M 2 . The sum of all possible re-scalings with subsequent expansions with non-trivial 
coefficients then reproduces the expansion of the result of the original integral. 
To be specific, let us apply this program to H: 



H(p 2 ,m 2 ,M 2 ;n) = -i / — — r = . 238 

J 2vr ) n k 2 M 2 ^ - 2p ■ kM a * + p 2 - m 2 + i0+ k 2 M 2 ^ — M 2 + i0+ 



From Eq. (12381) we see that the first fraction does not contribute to the overall factor M l3i for any 
It will be expanded in (positive) powers of (k 2 M 2ai — 2p ■ kM ai ) except for ctj = 0. For < on < 1, we 
re-write the second fraction as 

1 1 11/ M 2 ~ 2a > \ , . 

1 + T - — + • ■ ■ . (239) 



M 2oii (k 2 — M 2 ~ 2oti + z0+) M 2a >k 2 + i0+\ k 2 + i0 
On the other hand, if 1 < a, we re-write the second fraction as 

1 1 



— (l + k 2 M 2a *~ 2 + •••). (240) 



M 2 {k 2 M 2 ^- 2 -l + i0+) M 2 

In both cases one obtains integrals of the type / d n k k^ 1 ■ ■ ■ k^ m as the coefficients of the expansion. 
However, such integrals vanish in dimensional regularization. Therefore, the only non-trivial terms in 
the sum of Eq. (I237P correspond to either a,i = or = 1. Thus we obtain 

H(p 2 ,m 2 ,M 2 ;n) = H ( -°\p 2 , m 2 , M 2 ; n) + H < ' 1 \p 2 , m 2 , M 2 ; n), (241) 

where 

oo „ A n h. 1 1 

ff<°V, m\ M 2 ; n) = -ig («•)' / ^ - _ ^ - - ^ _ ^ - - ^ - ^ , (242, 

and 



M n_2 00 f-lV'M- 7 ' /■ ri n fr (k 2 M-2p-kY 
H^Mn) = -^E^/^^ (243) 

A comparison with the direct calculation of -ff shows that the dimensional counting method indeed leads 
to the correct expressions |Gegelia et al., 1994] . While the loop integrals of Eq. (12431) have a simple 



analytic structure in p 2 — m 2 , the same technique can be repeated for the loop integrals of Eq. (12421) 
when p 2 — m 2 — ► 0, now using the change of variable k = (p 2 — m 2 ) ll k with arbitrary non-negative real 
numbers ji. 



5 Applications 

In the following we will illustrate a few selected applications of the manifestly Lorentz-invariant frame- 
work to the one-nucleon sector. 
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Figure 17: Contributions to the nucleon self energy at 0(q 4 ). The number n in the interaction blobs 
refers to C^. The Lagrangian C^ N does not produce a contribution to the ttNN vertex. 

5.1 Nucleon mass and sigma term at 0(q 4 ) 

A full one-loop calculation of the nucleon mass also includes 0(q 4 ) terms (see Fig. [T71) . The quark-mass 
expansion up to and including 0(q 4 ) is given by 

m N = m + hM 2 + k 2 M 3 + k 3 M 4 In ( —] + k 4 M A + C(M 5 ), (244) 

\m J 

where the coefficients fcj in the EOMS scheme read Fuchs et al., 2003a| 



3g 2 3 
k\ = — 4ci, k 2 = k 3 = — - — ( g\ — 8cim + Coin + 4c 3 m ) , 

ki = 32^< 1+4Cim » + 128^ C2 -^- < 245 » 



Here, e\ = We^s + 2ens + 2en6 is a linear combination of C(g 4 ) coefficients |Fettes et al., 2 000 1. A 



comparison with the results using the infrared regularization |Becher and Leutwyler, 19 99 shows that 
the lowest-order correction (k\ term) and those terms which are non-analytic in the quark mass rh (k 2 
and k 3 terms) coincide. On the other hand, the analytic /c 4 term (~ M 4 ) is different. This is not 
surprising; although both renormalization schemes satisfy the power counting specified in Sec. 14.2.21 
the use of different renormalization conditions is compensated by different values of the renormalized 
parameters. 

For an estimate of the various contributions of Eq. (12441) to the nucleon mass, we make use of the 
parameter set 

c\ = — 0.9 m]^ 1 , C2 = 2.5m^r 1 , C3 = — 4.2 m^ 1 , C4 = 2.3m^ 1 , (246) 



which was obtained in Ref. Becher and Leutwyler, 2001 from a (tree-level) fit to the ttN scattering 



threshold parameters. Using the numerical values 

g A = 1-267, Ft, = 92.4 MeV, m N = m p = 938.3 MeV, = M w+ = 139.6 MeV, (247) 

one obtains for the mass of nucleon in the chiral limit (at fixed m s 7^ 0): 

m = m N - Am = [938.3 - 74.8 + 15.3 + 4.7 + 1.6 - 2.3 ± 4] MeV = (883 ± 4) MeV (248) 

with Am = (55.5 ± 4) MeV. Here, we have made use of an estimate for eiM 4 = (2.3 ±4) MeV obtained 
from the a term. (Note that errors due to higher-order corrections are not taken into account.) In 
terms of the SU(2)ixSU(2)/j-chiral-symmetry-breaking mass term of the QCD Hamiltonian, 

H sh = m(uu + dd) 1 (249) 

the pion-nucleon a term is defined as the proton matrix element 

^ = ^-(p(p,«)|Wsb(0)b(p J s)> (250) 



,- p 



66 



at zero momentum transfer. The a term provides a sensitive measure of explicit chiral symmetry 
breaking in QCD, because it is a correction to a null result in the chiral limit rather than a small 



correction to a non-trivial result Pagels, 1975 . The quark-mass expansion of the a term reads 

a x M 2 + a 2 M 3 + a 3 M 4 In ( —] + a A M A + C(M 5 ), (251) 

V m J 



a 



with 



a 4 



-4ci, 
3 



02 



9g, 



0"3 



8n 2 F 2 m 



3g 2 



64ttF 2 

2 

1 1 Cim(l + 2g A 



l6n 2 F 2 m 
2 , c 3 m 



kim + c 2 m 



+ 4c 3 m) 



2ei. 



We obtain [with e x = in Eq. fl252|) 



a = (74.8 - 22.9 - 9.4 - 2.0) MeV = 40.5 MeV. 



(252) 



(253) 



The result of Eq. (12531) has to be compared with, e.g., the dispersive analysis a = (45 ± 8) MeV of Ref. 
Gasser et al., 199~T] which would imply, neglecting higher-order terms, — 2eiM 4 « (4.5 ± 8) MeV. As 



has been discussed, e.g., in Ref. |Becher and Leutwyler, 1999] , a fully consistent description would also 
require to determine the low-energy coupling constant c\ from a complete 0(q A ) calculation of, say, ttN 
scattering. The results of Eqs. (I245[) and (12521) satisfy the constraints as implied by the application of 
the Hellmann-Feynman theorem to the nucleon mass |Gasser et al., 1988] , 



a 



iVF 



dm 



N 



dM 2 



(254) 



5.2 Chiral expansion of the nucleon mass to 0(q 6 ) 



So far, essentially all of the manifestly Lorentz-invariant calculations have been restricted to the 
one-loop level. One of the exceptions is the chiral expansion of the nucleon mass which, in the 
framework of the reformulated infrared regularization, has been calculated up to and including 0(q 6 ) 



Schindler et al., 2007b, Schindler et al., 2008 



m N = m + hM 2 + k 2 M 3 + k 3 M 4 In — + /,••.!/ : + k 5 M 5 In — + k 6 M 5 + k 7 M 6 In 2 — + k 8 M 6 In — + k 9 M 6 . 

fX fj, II fl 

(255) 

We refrain from displaying the lengthy expressions for the coefficients fcj but rather want to discuss 
a few general implications [Schindler et al.,"20 08|. Chiral expansions like Eq. (12551) currently play an 
important role in the extrapolation of lattice QCD results to physical quark masses. Unfortunately, 
the numerical contributions from higher-order terms cannot be calculated so far since, starting with 
&4, most expressions in Eq. (12551) contain unknown low-energy coupling constants (LECs) from the 
Lagrangians of 0(q 4 ) and higher. The coefficient kr, is free of higher-order LECs and is given in terms 
of the axial-vector coupling constant g A and the pion-decay constant F: 



fen 



3g? 



1024tt 3 F 4 



(leg 2 ! - s) . 



While the values for both g A and F should be taken in the chiral limit, we evaluate k% using the 
physical values gA = 1.2695(29) and F n = 92.42(26) MeV. Setting (j, = mjv, = {rn p + m n )/2 = 
938.92 MeV, and M = = 139.57 MeV we obtain k 5 M 5 \n(M/m N ) = -4.8 MeV. This amounts 
to approximately 31% of the leading non-analytic contribution at one- loop order, fc 2 M 3 . Figure [18] 
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Figure 18: Pion mass dependence of the term k§M 5 ln(M/mjv) (solid line) for M < 400 MeV. For 
comparison also the term k 2 M 3 (dashed line) is shown. 



shows the pion mass dependence of the term k 5 M 5 ln(M/mjv) (solid line) in comparison with the 
term /c 2 M 3 (dashed line) for pion masses below 400 MeV which is considered a region where chiral 
extrapolations are valid (see, e.g., Refs. |Meifiner, 2006| , |Djukanovic et al., 2006| ). We see that already 
at M « 360 MeV the term /c 5 M 5 h\(M/m^) becomes as large as the leading non-analytic term at one- 
loop order, k 2 M 3 , indicating the importance of the fifth-order terms at unphysical pion masses. Our 
results for the renormalization-scheme-independent terms agree with the heavy-baryon ChPT results of 



Ref. McGovern and Birse, 1999 



5.3 Form factors of the nucleon 
5.3.1 Scalar form factor 

The pion-nucleon a term corresponds to the kinematical point t = of the scalar form factor which is 
defined as 

W,s')\H sh (0)\p(p,s)) =u( P ',s')u(p,s)a(t), t={p'-p) 2 . 

The numerical results for the real and imaginary parts of the scalar form factor at 0(q 4 ) are shown 
in Fig. [19] for the extended on-mass-shell scheme (solid lines) and the infrared regularization scheme 
(dashed lines). While the imaginary parts are identical in both schemes, the differences in the real parts 
are practically indistinguishable. Note that for both calculations <r(0) and A CT = a(2M%) — cr(0) have 
been adjusted to the dispersion results of Ref. Gasser et al., 1991] , A CT = (15.2 ± 0.4) MeV. 

Figure [20] contains an enlargement near t ~ 4M^ for the results at 0(q 3 ) which clearly displays how 
the heavy-baryon calculation fails to produce the correct analytic behavior not only at the tree level 
but also in higher-order loop diagrams. Both real and imaginary parts diverge as t — > 4M%. 



5.3.2 Electromagnetic form factors 

Imposing the relevant symmetries such as translational invariance, Lorentz covariance, the discrete sym- 
metries, and current conservation, the nucleon matrix element of the electromagnetic current operator 

2 1 - 

JP(x) = -u(x)^u(x) - -d(x)^d(x), 
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Figure 19: Scalar form factor a(t) as a function 
of t at C(g 4 ). 



Figure 20: Real and imaginary parts of the 
scalar form factor CIS db function of t at 0(q 3 ) 
in the vicinity of t = AM 2 . Solid lines: EOMS 
scheme; dashed lines: infrared regulariza- 
tion (IR) of Ref. Becher and Leutwyler, 19 99 1; 
dotted lines: HBChPT calculation of Ref. 
|Bernard et al., 1992a . On this scale the (un- 
physical) divergence of both real and imagi- 
nary parts of the heavy-baryon result becomes 
visible. 



can be parameterized in terms of two form factors, 



{N{p',s')\J^0)\N(p,s))=u(p',s') 



F l N (Q 2 )^ + i^F 2 N (Q 2 ) 



u(p,s), N = p,n, 



(256) 



where q = p' — p, Q 2 = —q 2 , and m p is the proton mass. At Q 2 = 0, the so-called Dirac and Pauli form 
factors Fi and F 2 reduce to the charge and anomalous magnetic moment in units of the elementary 
charge e and the nuclear magneton e/(2m p ), respectively, 

Ff (0) = 1, F™(0) = 0, F$(0) = 1.793, F 2 n (0) = -1.913. 

The Sachs form factors Ge and Gm are linear combinations of F\ and F 2 , 



g e(Q ) = F"(Q 



4m? 



(Q*), GZ{Q 2 )=F[\Q 2 ) + Fi\Q 2 ), N = p,n 



and, in the non-relativistic limit, their Fourier transforms are commonly interpreted as the distribution of 
charge and magnetization inside the nucleon. For a covariant interpretation in terms of the transverse 
charge density see Refs. Miller, 2007 , Carlson and Vanderhaeghen, 2008 . The description of the 



electromagnetic form factors of the nucleon presents a stringent test for any theory or model of the 
strong interactions (see, e.g., Ref. |Perdrisat et al., 2007 for a recent review). 

In the framework of chiral perturbation theory, the electromagnetic form factors were calculated in 
the early relativistic approach |Gasser et al., 1988] , the heavy-baryon approach |Bernard et al., 199 2a|, 
Fearing et al., 1997| , the small-scale expansion Bernard et al., 19 98 1 , the infrared regularization 



Kubis and MeiBner, 2001] , and the EOMS scheme Fuchs et al., 2004 . All these calculations have in 



common that they fail to describe the proton and nucleon form factors for momentum transfers beyond 
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Figure 21: The Sachs form factors of the nucleon in manifestly Lorentz-invariant chiral pertur- 
bation theory at 0(q A ) without vector mesons. Full lines: results in the extended on-mass-shell 
scheme; dashed lines: results in infrared regularization. The experimental data are taken from Ref. 
Friedrich and Walcher, 2003| . 



Q 2 ~ 0.1 GeV 2 . Moreover, up to and including 0(q 4 ), the most general effective Lagrangian provides 
sufficiently many independent parameters such that the empirical values of the anomalous magnetic 
moments and the charge and magnetic radii are fitted rather than predicted. Figure [211 shows the Sachs 
form factors in the momentum transfer region GeV 2 < Q 2 < 0.4 GeV 2 in the EOMS scheme and the 



reformulated infrared regularization Schindler et al., 2005| . 
In Ref. " 



Kubis and Meifiner, 2001 



it was shown that the inclusion of vector mesons can result in the 
re-summation of important higher-order contributions. In standard ChPT, such vector meson contri- 
butions manifest themselves in terms of the values of the low-energy coupling constants. Symbolically, 
the contributions to certain LECs originate from the expansion of the vector-meson propagator, 



1 



1 



q 2 -M 2 



M 2 



M 2 



Ml, 



0(q 6 



combined with the relevant vector-meson vertices. However, diagrams with internal vector-meson lines 
inside loops were not considered, because a generalization of ChPT which fully includes the effects of 
vector mesons as intermediate states in loops was not yet available Kubis and Meifiner, 2001 . On the 
other hand, the EOMS renormalization scheme of Ref. IFuchs et al., 2003a| and the reformulated version 



of infrared regularization of Ref. |Schindler et al., 2004a| both allow to include virtual vector mesons 



systematically in the region of the applicability of baryon chiral perturbation theory |Fuchs et al., 20 03b 



(see also Ref. |Bruns and Meifiner, 2008] ). The standard power counting determines which diagrams 
(including diagrams with vector mesons appearing in loops) should be taken into account to a given 
order in the chiral expansion. 

In Ref. [ Schindler et al., 2005| the electromagnetic form factors were calculated with the p, u, and 
mesons as explicit degrees of freedom. In the vector-field representation of Ref. |Ecker et al., 1989b] the 
p meson is represented by p^ = p^Ti and the u and mesons by uj^ and M , respectively. The coupling 
of the vector mesons to pions and external fields is at least of 0(q 3 ), 



A3) 



where the field strength tensors are given by 



fM^O - f^ u f! s) 



I at 



+ 



(257) 
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(I) (II) 




(V) 

Figure 22: Feynman diagrams including vector mesons that contribute to the electromagnetic form 
factors of the nucleon up to and including 0(q A ). External leg corrections are not shown. Solid, wiggly, 
and double lines refer to nucleons, photons, and vector mesons, respectively. The numbers in the 
interaction blobs denote the order of the Lagrangian from which they are obtained. The direct coupling 
of the photon to the nucleon is obtained from C^ N and C^n- 



with f^ u and defined in Eqs. (11 111) and (11121) . respectively. For the case of a coupling to an external 
electromagnetic potential A^, the external fields are given by Eq. (I5TI) . Furthermore, in terms of the 
connection V 1 of Eq. (11781) . we define 

p ^ = vy - w, vv = <9V + [r M , p v \ , 

and, finally, 

oj^ = d»u v - <9V, 0^ = &*<f> v - d v ^. 
The lowest-order Lagrangian for the coupling to the nucleon is given by 

C vn = \ E gvV-fVfl, (258) 

V=p,LO,rj> 

and the 0(q) Lagrangian reads 

C vn = \ E Gv^a^V^. (259) 

V=p,LU,(j) 

The additional power-counting rules state that vertices from C^y count as 0(q 3 ) and vertices from dy N 
as 0(q l ), respectively, while the vector-meson propagators count as O(q ). The additional diagrams 
involving vector mesons that contribute in the calculation of the form factors up to and including 
0(q 4 ) using the Lagrangians of Eqs. (125 7p . (12581) . and (I259j) are shown in Fig. [22J The parameters 
of the vector-meson Lagrangian of Eq. (12571) for the coupling to external fields have been taken from 
Ref. |Ecker et al., 1989b] , and those of Eqs. (I258P and (12591) for the coupling of vector mesons to the 



nucleon from the dispersion relations of Refs. |Mergell et al., 1996 , |Hammer and Meifiner, 2004 . 
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Figure 23: The Sachs form factors of the nucleon in manifestly Lorentz-invariant chiral perturbation 
theory at 0(q 4 ) including vector mesons as explicit degrees of freedom. Full lines: results in the extended 
on-mass-shell scheme; dashed lines: results in infrared regularization. The experimental data are taken 



from Ref. Friedrich and Walcher, 2003 . 



As expected on phenomenological grounds, the quantitative description of the data has improved 
considerably for Q 2 > 0.1 GeV 2 (see Fig. 123]) . The small difference between the two renormalization 
schemes is due to the way how the regular higher-order terms of loop integrals are treated. Numer- 
ically, the results are similar to those of Ref. |Kubis and Meifiner, 2001] . Due to the renormalization 
condition, the contribution of the vector-meson loop diagrams either vanishes (IR) or turns out to be 
small (EOMS). Thus, in hindsight our approach puts the traditional phenomenological vector-meson- 
dominance model on a more solid theoretical basis. In the sense of a strict chiral expansion in terms of 
small external momenta q and quark masses m q at a fixed ratio m q /q 2 |Gasser and Leutwyler, 1984] , up 
to and including C(g 4 ) the results with and without explicit vector mesons are completely equivalent. 
The additional vector- meson contributions up to this order are compensated by a readjustment of the 
low-energy constants pertaining to the theory including vector mesons as dynamical degrees of freedom. 
On the other hand, the inclusion of vector-meson degrees of freedom in the present framework results 
in a reordering of terms which, in an ordinary chiral expansion, would show up at higher orders beyond 
0{q i ). It is these terms which change the form factor results favorably for larger values of Q 2 . It 
should be noted, however, that this re-organization proceeds according to well-defined rules so that a 
controlled, order-by-order, calculation of corrections is made possible. In contrast to the calculation 
without vector mesons, the Sachs form factors G P E , G P M , and now show sufficient curvature to gen- 
erate a more accurate phenomenology for values of Q 2 , where the ordinary chiral expansion to the same 
order is no longer reliable. 



5.3.3 Axial and induced pseudoscalar form factors 

Assuming isospin symmetry, the most general parametrization of the isovector axial-vector current 
evaluated between one-nucleon states is given by 



(N(p',s')\A?(0)\N(p,s))=u(p',s') 



Yl,G A {Q 2 ) + -^ 75 G P (g 2 ) 
2m N 



^u(p,s), (260) 



where q = p' — p, Q 2 = —q 2 , and mjy denotes the nucleon mass. Ga{Q 2 ) is called the axial form 
factor and Gp(Q 2 ) is the induced pseudoscalar form factor. The value of the axial form factor at zero 
momentum transfer is defined as the axial-vector coupling constant, qa = Ga(Q 2 = 0) = 1.2695(29) 
Amsler et al., 2008] , and is quite precisely determined from neutron beta decay. The Q 2 dependence 
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of the axial form factor can be obtained either through neutrino scattering or pion electroproduction 
Bernard et al., 2002 1 and references therein). The second method makes use of the so-called 



sec 



Adler-Gilman relation Adler and Gilman, 1966] which provides a chiral Ward identity establishing a 



connection between charged pion electroproduction at threshold and the isovector axial-vector current 



evaluated between single-nucleon states (see, e.g., Ref. |Scherer and Koch, 1991| for more details). The 
induced pseudoscalar form factor Gp{Q 2 ) has been investigated in ordinary and radiative muon capture 
as well as pion electroproduction (see Ref. |Gorringe and Fearing, 2004| for a review). 

For the analysis of experimental data, Ga{Q 2 ) is conventionally parameterized using a dipole form 

as 



Ga(Q 2 ) = - 9A n2 „ , (261) 



A 



where the axial mass M A is related to the axial root-mean-square radius by (r\) 2 = 2y/%/M A . The global 
average for the axial mass extracted from neutrino scattering experiments given in Refs. 
|Liesenfeld et al., 1999| , |Bernard et al., 2002| is 



M A = (1.026 ± 0.021) GeV. (262) 

The extraction of the axial mean-square radius from charged pion electroproduction at threshold is 
motivated by the current algebra results and the PCAC hypothesis. At threshold (the spatial compo- 
nents of) the center-of-mass transition current for pion electroproduction can be written in terms of 
two s-wave amplitudes E 0+ and L Q+ , 



eM 



thr rriN 



4ttW r 

ia±E +(k 2 ) + ia\\L 0+ (k 2 



where W is the total center-of-mass energy, k 2 is the four momentum transfer squared of the virtual 
photon, and a\\ = a ■ kk and a± = a — B\\. The reaction p(e, e'ir + )n has been measured at MAMI at an 
invariant mass of W = 1125 MeV (corresponding to a pion center-of-mass momentum of |<f | = 112 MeV) 
and photon four-momentum transfers of —k 2 = 0.117,0.195 and 0.273 GeV 2 Liesenfeld et al., 19 99 1. 
Using an effective-Lagrangian model an axial mass of 

M A = (1.077 ±0.039) GeV 

was extracted, where the bar is used to distinguish the result from the neutrino scattering value. In the 
meantime, the experiment has been repeated including an additional value of — k 2 = 0.058 GeV 2 and is 
currently being analyzed. The global average from several pion electroproduction experiments is given 
by |Bernard et al., 2002| 

M A = (1.069 ± 0.016) GeV. (263) 

It can be seen that the values of Eqs. f 1 2 6 2 j) for the neutrino scattering experiments are smaller than 
Eq. (12631) for the pion electroproduction experiments. The discrepancy was explained in heavy-baryon 
chiral perturbation theory |Bernard et al., 1992b| . It was shown that at 0(q 3 ) pion loop contribu- 
tions modify the k 2 dependence of the electric dipole amplitude from which M A is extracted. These 
contributions result in a change of 

AM A = 0.056 GeV, (264) 

bringing the neutrino scattering and pion electroproduction results for the axial mass into agreement. 
In a recent analysis |Bodek et al., 2008] updated expressions for the vector form factors have been 
taken into account together with the hadronic correction of Eq. (12641) to produce an average from both 
neutrino and electroproduction experiments, 

M A = (1.014 ± 0.014) GeV. (265) 
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Figure 24: The axial form factor Ga in manifestly Lorentz-invariant ChPT at 0(q A ). Full line: result in 
infrared renormalization with parameters fitted to reproduce the axial mean-square radius corresponding 
to the dipole parametrization with Ma = 1.026 GeV (dashed line). The dotted and dashed-dotted 
lines refer to dipole parameterizations with Ma = 0.95 GeV and Ma = 1.20 GeV, respectively. The 
experimental values are taken from Bernard et al., 2002 1. 



Earlier calculations of the axial form factor were performed in the framework of heavy-baryon ChPT 
Bernard et al., 1992b , Fearing et al., 1997 and the small-scale expansion Bernard et al., 1998|. In 



Ref. Schindler et al., 2007a the form factors Ga and Gp have been calculated in manifestly Lorentz- 



invariant baryon ChPT up to and including 0(q 4 ). The axial form factor can be written as 

G A (Q 2 ) =9A-\g A (ri> Q 2 + ^r 2 L(Q 2 ), (266) 

where (r A ) is the axial mean-square radius and L contains loop contributions and satisfies L(0) = 
1/(0) = 0. The result for Ga in the momentum-transfer region GeV 2 < Q 2 < 0.4 GeV 2 is shown 
in Fig. EH The parameters have been determined such as to reproduce the axial mean-square radius 
corresponding to the dipole parametrization with Ma = 1.026 GeV (dashed line). The dotted and 
dashed-dotted lines refer to dipole parameterizations with Ma = 0.95 GeV and Ma = 1-20 GeV, 
respectively. The loop contributions from L(Q 2 ) are small and the result does not produce enough 
curvature to describe the data for momentum transfers Q 2 > 0.1 GeV 2 . The situation is similar to the 
electromagnetic case of Fig. EH where ChPT at 0(q A ) also fails to describe the form factors beyond 
Q 2 > 0.1 GeV 2 . 

In addition to the standard treatment including the nucleon and pions, the axial-vector meson 
ai(1260) has also been considered as an explicit degree of freedom |Schindler et al., 2007a . In the 



vector-field formulation of |E cker et al., 1989b] the ai(1260) meson is represented by = A^ti. The 
advantage of this formulation is that the coupling of the axial-vector mesons to pions and external 
sources is at least of 0(q 3 ). The calculation of the contributions to the isovector axial- vector form 
factors only requires the term 

4 3 ] = ^TrOV/r), (267) 
where the field strength tensor is defined as 

r = u + f%u-uf^u\ 
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Figure 25: Diagram containing the axial- vector meson (double line) contributing to the form factors at 



and 

Apt, = V fj,A u — T \/ u A fM , V ^Ay = d^A v + [T^, A u ], 

with the connection of Eq. (j!78p . The coupling of the axial- vector meson to the nucleon starts at O(q ). 
The corresponding Lagrangian reads 

4°i = ^YlsA^. (268) 

A calculation up to and including 0(q 4 ) would in principle also require the Lagrangian of 0(q). However, 
there is no term at this order that is allowed by the symmetries. The additional power-counting rules 
are as in Section [5.3.21 for the vector mesons. We count the axial- vector meson propagator as O(q ), 
vertices from as 0(q 3 ), and vertices from Can as O(q ), respectively. The contributions of the 
axial- vector meson to the form factors G a and Gp at 0(q A ) originate from the diagram in Fig. [23 The 
inclusion of the axial-vector meson effectively results in one additional low-energy coupling constant 
which has been determined by a fit to the data for Ga{Q 2 )- The inclusion of the axial- vector meson 
results in an improved description of the experimental data for Ga (see Fig. 126]) . while the contribution 
to Gp is small. 



5.3.4 Pion-nucleon form factor 

The pion-nucleon form factor G n N{Q 2 ) may be defined in terms of the pseudoscalar quark density 



P% — iqibftq and the average light-quark mass m as [Gasser et al., 1988 



A/f 2 P 

m(N( P ',s')\Pm\N(p,s)) = M / + q 2 GMQ 2 )W, aO Wfo s l ( 269 ) 

where q = p' — p, Q 2 = — q 2 , and $j(x) = ^|f) j s the corresponding interpolating pion field. The 
pion-nucleon coupling constant is given by g^N = G n N(—M%). Using the (QCD-) partially conserved 
axial- vector current (PCAC) relation, d^Af = rhPi, the pion-nucleon form factor is completely given in 
terms of the axial and the induced pseudoscalar form factors, 

This is an exact relation which holds true for any value of Q 2 . The result at 0(q 4 ) is given by 



Schindler et al., 2007a 



n m N9A a Q 2 

G nN (Q ) = — 9,n^j^ + ■■■ 

where A = 1 — ™ N g 9 f r denotes the Goldberger-Treiman discrepancy. The chiral expansion of the pion- 



nucleon coupling constant can be found in Ref. Schindler et al., 2007a 
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Figure 26: Left panel: Axial form factor Ga in manifestly Lorentz-invariant ChPT at 0(q 4 ) including 
the axial- vector meson ai(1260) explicitly. Full line: result in infrared renormalization, dashed line: 
dipole parametrization. The experimental data are taken from Ref. Bernard et al., 2002] . Right panel: 
The induced pseudoscalar form factor Gp in manifestly Lorentz-invariant ChPT at 0(q 4 ) including the 
axial- vector meson ai(1260) explicitly. Full line: result with axial- vector meson; dashed line: result 
without axial- vector meson. One can clearly see the dominant pion pole contribution at Q 2 ~ —M 2 . 



6 Conclusion 

Effective field theory has become a very important tool for investigating the dynamics of the strong inter- 
actions. In particular, mesonic chiral perturbation theory is a full-grown and mature area of low-energy 
particle physics which has successfully been applied at the two-loop level. Whether the predictions for 
the electromagnetic polarizabilities of the charged pion are really in conflict with empirical data re- 
mains to be seen. In the baryonic sector new renormalization conditions have reconciled the manifestly 
Lorentz-invariant approach with the standard power counting. Phenomenological extensions allowing 
for the rigorous inclusion of (axial-) vector-meson degrees of freedom (and also of the A(1232) resonance) 
have opened the door to an extended kinematic region. Unfortunately, the question of convergence in 
the three-flavor sector remains a controversial issue |Lehnhart et al., 2005 , even though the manifestly 
Lorentz-invariant approach might yield better phenomenological results Geng et al., 2008] . Finally, be- 



yond the one-nucleon sector the covariant framework has been used in the discussion of relativistic cor- 
rections to the nucleon-nucleon potential (see, e.g., Refs. [Higa and Robilotta, 2003| , |Robilotta, 2007| ) 



or may be applied to the nuclear many-body problem (see, e.g., Refs. |Furnstahl, 2004| , |Serot, 2004| ). 
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